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Abstract 


INTERATOMIC  FORCES  IN  TRANSITION  METALS 

In  this  work  we  have  been  concerned  with  obtaining  a 
practical  description  of  the  force  fields  in  transition 
metals.  There  is  some  evidence  that  directional  effects 
due  to  the  d  orbitals  of  the  transition  metals  increase 
through  the  sequence  face  centred  cubic  (f.c.c.),  hexagonal 
close  packed  (h.c.p.)  and  bofly  centred  cubic  (b.c.c.) 
structures.  '  These  features  have  led  us  to  study  the  utility 
of  density-independent  central  pair  forces,  plus  a  volume- 
dependent  energy  for  the  f.c.c.  metals  Ni  and  Pt.  We  have 
constructed  central  pair  potentials  for  both  these  metals 
which  give  fair  agreement  between  the  calculated  and  exper¬ 
imentally  determined  phonon  dispersion  curves  and  liquid 
structure  factors.  While  the  above  results  strongly  demon¬ 
strate  the  usefulness  of  such  a  representation  of  the  force 
fields  in  the  f.c.c.  solid  transition  metals,  we  have,  by 
studying  the  X-ray  scattering  factors,  found  the  charge  dist¬ 
ribution  in  the  b.c.c.  transition  metals  to  be  highly  angular. 
Using  a  Wannier  function  formalism,  we  have  found  that  this 
angularity  can  be  well  represented  by  a  bond  charge  model  in 
which  the  bond  charges  are  situated  at  nearest  and  next  nearest 
neighbour  interstitials.  The  angularity  in  the  charge  dist¬ 
ribution  is  a  result  of  degenerate  bands  forming  hybrid  wave- 
functions.  In  h.c.p.  Be  the  degenerate  bands  lead  to  sp 
hybridization.  The  usefulness  of  the  Wannier  formalism  was 
investigated  quantitatively  by  calculating  the  X-ray  scatt¬ 
ering  factors  of  Be.  Expressions  for  the  elastic  constants 
and  phonon  dispersion  relations  have  been  derived  for  this 
bond  charge  model  and  used  to  calculate  these  properties  of 
Fe.  The  fair  agreement  with  the  experimentally  determined 
phonon  dispersion  curves  shows  that  in  addition  to  the  central 
ion-ion  pair  interactions,  there  also  exists  bond  electron-ion 
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CHAPTER  1 


Introduction 


Any  real  understanding  of  the  physics  of  metals  must 
entail  some  knowledge  of  the  force  fields.  Since  the  force 
fields  are  determined  by  the  interaction  of  the  ions  and 
electrons  amongst  themselves,  an  accurate  knowledge  of  the 
force  fields  can  be  obtained  if  the  electronic  distribution 
is  known.  Thus,  we  may  determine  the  force  fields  either 
directly,  by  analysing  experimental  data  (e.g.  the  phonon 
dispersion  data)  or  indirectly,  by  building  a  picture  of 
the  charge  density  in  the  metal. 

The  pseudopotential  theory  of  metals  which  makes  use 
of  the  free  electron  model  allows  one  to  write  the  total 
energy  of  a  crystal  in  terms  of  the  dielectric  function, 
which  describes  the  electronic  distribution.  For  a  free 
electron  gas,  and  working  to  second  order  in  the  perturbing 
electron- ion  interaction,  it  turns  out  that  this  may  be 
further  simplified  into  a  sum  of  effective  ion-ion  central 
pair  interactions  and  a  volume  dependent  structure  indepen¬ 
dent  energy.  Since  this  theory  has  proved  quite  successful 
for  the  simple  metals,  it  implies  that  the  force  fields  in 
these  metals  can  be  represented  by  central  pair  ion-ion 
interactions . 

Harrison  (1969)  extended  this  model  for  the  d-band 
metals  by  including  in  the  total  wavefunction  d-resonance 
states  to  account  for  the  unfilled  d-band,  and  also  a 
hybridization  term.  Using  this  model,  Moriarty  (1972)  was 
able  to  determine  effective  interionic  potentials  in  the 
noble  metals  by  working  to  second  order  in  the  perturbation. 
Although  hybridization  was  included,  the  pair  interactions 
were  taken  to  be  central.  This  approach  is  inconsistent 
with  the  angularity  present  in  the  charge  density  due  to 
hybridization. 


Bertoni  et  al.  (1975)  calculated  the  phonon  dispersion 
curves  and  elastic  constants  of  Be  and  found  that  an  exten- 
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sion  of  the  pseudopotential  theory  to  third  order  was 
needed  to  obtain  reasonable  agreement  with  experiment. 

Since  this  extension  yields  not  only  pair  forces,  but  aLso 
three  body  forces,  it  is  indicative  of  the  non-central 
nature  of  the  force  fields.  It  is  interesting  to  note  that 
the  main  feature  common  to  Be  and  the  d-band  metals  is  that 
of  overlapping  valence  bands.  It,  therefore,  appears  that; 
unlike  in  the  simple  metals,  the  force  fields  in  crystals 
with  degenerate  bands  are  non-central.  In  this  work,  we 
have  addressed  ourselves  to  two  questions: 

(i)  Although  the  force  fields  in  the  d-band  metals 
will  undoubtedly  be  non-central,  can  the  approx¬ 
imation  of  central  forces  in  some  d-band  metals 
still  be  a  useful  one. 

(ii)  In  these  d-band  metals  where  the  central  approx¬ 
imation  fails,  how  can  the  non-central  character 
of  the  forces  be  usefully  represented. 

A  first  principles  calculation  of  the  force  field  in  a 
metal  requires  knowledge  of  the  dielectric  function  of  the 
electron  distribution.  As  even  in  the  simple  metals  the 
uncertainty  in  the  dielectric  function  leads  to  widely 
varying  interatomic  potentials  (Ashcroft  and  Langreth  1967, 
Ho  1970)  there  is  little  hope  of  obtaining  a  good  potential 
for  any  of  the  transition  metals.  We  have  adopted,  there¬ 
fore,  a  procedure  which  is  semiempirical,  but  in  which  the 
essential  guidelines  laid  down  by  the  electron  theory  of 
d-band  metals  are  taken  into  account  as  follows: 

(i)  The  intermediate  and  long  range  interactions  are 

oscillatory  in  nature  and  have  the  form  cos(ar+b)/r 

(ii)  The  d-shells  lead  to  considerably  harder  repulsive 
core-core  interactions  than  in  the  simple  sp  metals 

(iii)  The  directional  effects  due  to  the  d-electrons 

appear  to  vary  according  to  the  structure  of  the 
crystalline  transition  metals.  If  one  adopts  an 
approach  based  on  hybrids  of  d  orbitals  with  s  and 
p  orbitals,  then  the  crystal  symmetry  indicates 
which  are  the  appropriate  hybrids. 
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These  guidelines  have  been  taken  into  account  in  our 
work,  in  an  admittedly  incomplete  way,  as  follows: 

(a)  It  is  a  well-known  fact  that  the  volume  change 
of  a  metal  upon  melting  is  sufficiently  small  to 
allow  one  to  extract  information  on  the  interatomic 
forces  in  the  solid  from  the  liquid  state.  In  a 
classical  liquid,  the  structure  factor  is  directly 
related  to  the  effective  interaction  potential. 

As  the  measured  liquid  structure  factors  of  the 
d-band  metals  do  not  differ  qualitatively  from 
those  of  the  simple  metals  (Huijben  and  van  der 
Lugt  1979,  Waseda  and  Ohtani  1974),  we  have 
assumed  that  in  the  liquid  state  the  force  fields 
are  also  qualitatively  the  same. 

(b)  From  liquid  metal  data  interpreted  on  the  basis  of 
pair  interactions ,  we  have  been  able  to  extrapolate 
approximate  interatomic  pair  potentials  at  near 
equilibrium  distances. 

(c)  Both  liquid  structure  and  radiation  damage  data 
have  been  used  to  obtain  an  estimate  of  core  hard- 

•  ness.  Also,  the  melting  curves  of  the  d-band 

metals  under  pressure  lends  support  to  our  estim¬ 
ates  of  core  repulsion. 

(d)  From  the  crystal  structure,  using  group  theoretical 
techniques,  appropriate  hybrids  have  been  constructed 
and  the  resulting  charge  distribution  modelled  by 
bond  charges. 

Having  obtained  suitable  representations  for  the  charge 
distribution  and  force  fields  in  particular  d-band  metals, 
we  have  proceeded  to  test  their  usefulness  by  examining 
their  effects  on  the  second  and  third-order  elastic  const¬ 
ants  by  studying  the  phonon  dispersion  curves  and  by  invest¬ 
igating  the  single  crystal  X-ray  scattering  experiments. 


In  Chapter  2  we,  by  representing  the  force  fields  in 
the  face  centred  cubic  (f.c.c.)  d-band  metals  by  central 
pair  interactions,  calculate  pair  potentials  for  Ni 
and  Pt  using  liquid  structure  and  phonon  dispersion  data 
(Matthai  et  al.  1978a). 


The  angularity  in  the  charge  distribution  in  the 
transition  metals  is  discussed  in  Chapter  3,  and  a  form¬ 
alism  for  the  X-ray  scattering  factors  in  solids,  in  terms 
of  Wannier  functions  is  set  up  (Matthai  et  al.  1978b). 

With  this  formalism,  the  measurements  of  the  body  centred 
cubic  (b.c.c.)  transition  metals  are  explained.  Since  a 
quantitative  explanation  requires  us  to  know  the  Wannier 
functions,  we,  by  calculating  them  approximately,  have 
determined  the  X-ray  scattering  factors  of  Be.  This  is 
presented  in  Chapter  4. 

In  Chapter  5,  by  means  of  a  model  the  angularity  of 
the  charge  distribution  is  incorporated  into  the  force 
fields.  This  enables  us  to  calculate  the  phonon  dispersion 
curves  of  Fe.  Finally,  the  accuracy  of  the  potentials 
presented  in  this  work  are  discussed,  and  ways  to  improve 
them  are' suggested  in  Chapter  6 . 


CHAPTER  2 


Pair  potentials  in  the  f.c.c.  transition  metals 
§2. 1  Introduction 


If  we  consider  a  metal  to  be  a  collection  of  ions  pos¬ 
itioned  at  lattice  sites  surrounded  by  some  inhomogeneous 
electron  distribution,  then,  we  can  write  the  total  energy, 
E,  of  the  metal  as 

E  =  Ec  +  EA  +  Ev  (2.1) 

where  E^  +  Ey  is  that  part  arising  from  only  considering 
the  spatially  homogeneous  part  of  the  electron  distribution 
and  E  I*  is  the  extra  energy  arising  from  the  spatial  inhomo- 
geniety.  From  experimental  evidences  about  the  charge 
density  distributions,  it  appears  that  the  charge  distrib¬ 
utions  of  the  f.c.c.  transition  metals  are  more  homogeneous 
than  that  of  the  b.c.c.  metals.  We  have  therefore  made 
the  approximation  that  E^  may  be  neglected  for  the  f.c.c. 
metals  and  thus  can  be  treated  like  a  simple  metal  (sp). 

In  the  simple  metals  the  model  can  be  simplified  to 
one  of  spherically  symmetric  charged  ions  situated  at 
lattice  sites  immersed  in  a  homogeneous  electron  gas.  Thus 
the  ions  interact  between  themselves  directly  and  indirectly 
via  the  electron  gas.  The  contribution  to  the  energy  from 
the  direct  ion-ion  interaction  consists  of  the  direct  coul¬ 
omb  repulsion  between  the  ions,  and  the  energy  arising  from 
the  repulsion  due  to  the  overlap  of  the  ion  cores,  Ep.  In 
the  theory  of  the  simple  metals,  it  can  be  shown  (Harrison 
1966)  that  the  direct  coulomb  repulsion  together  with  the 
screened  interaction  gives  rise  to  a  contribution  to  the 
total  energy  of  Ep  +  Ey ,  where 

Ep  =  [  ♦<r..)  (2.2); 

*  lj  13 

$(r)  =  — 3  cos  2kpr  (2.3), 

Ey  is  a  volume  dependent  structure  independent  energy  and 
kp  is  the  Fermi  wavevector. 


The  change  in  energy  in  bringing  two  ions  from  infinity 
to  a  separation  R,AE,  can  be  shown  to  be  (Dick  and 
Overhauser  1958) 

AE  =  eqex  ;  qgx  =  4S2  (2.4) 

~2R" 

where  S  is  the  overlap  integral.  Since  this  energy  is 
positive,  there  is  overall  repulsion  between  the  atoms. 

Over  the  range  under  consideration,  this  change  in  energy 
can  also  be  fitted  by  the  Born-Mayer  (B-M)  form,  exp(-R/RQ) 
and  for  more  complex  shelled  ions  we  expect  the  same  form 
for  the  repulsive  energy,  namely 

AE  =  Y<lex  =  A  exp(-bR)  (2.5) 

R 

where  y,  A  and  b  are  constants. 

In  the  rare  gas  solids,  the  core-core  repulsion  is 
often  characterised  by  an  inverse  power  potential,  as  in 
the  case  of  the  Lennard-Jones  potential. 

The  repulsive  potential  expressed  either  in  the  B-M 
or  inverse  power  form 

<j>R(r)  =  A  exp(-br)  =  e(^)n  (2.6) 

with  A,  b,  e  and  o  as  parameters,  gives  the  contribution 
Ep  to  the  total  energy: 

E.=i  [  <2.7) 

R  R.*0  R  D 

Both  the  forms  given  above  have  been  used  to  explain 
experimental  results  (Gibson  et  al.  1960,  Lennard-Jones  1924). 
To  correlate  these  two  forms  and  for  some  of  the  transition 
metals,  we  have  examined  the  radiation  damage  data  and 
the  melting  curves  of  these  metals. 

Thus  assuming  the  inverse  power  form  for  E^,  the  total 
energy,  in  the  simple  metals  takes  the  form 
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E 

sm 


EC  +  EV 


=  l  *eff(rij)  +  EV(V)  (2‘8> 

with 

$eff-(r)  =  —  +  cos  (2kfr  +  6)  (2.9) 

r  r 

For  the  f.c.c.  metals  we  have  assumed  the  energy  to 
be  given  by  equation  (2.8)  with  the  interionic  potential 
having  the  same  form  as  in  equation  (2.9).  We  have,  how¬ 
ever,  treated  the  constants  in  the  potential  as  parameters. 

Since  we  have  assumed  the  ions  to  be  spherically 
symmetric,  upon  melting,  the  metals  acts  like  a  classical 
liquid.  Also  since  we  do  not  expect  the  electronic  dist¬ 
ribution  to  change  a  great  deal,  the  particles  of  the 
liquid  would  interact  with  each  other  through  a  potential 
that  would  be  very  nearly  equal  to  However,  the 

disorder  present  in  the  liquid  would  result  in  a  damping 
of  pair  potential  by  a  factor  exp(-Xr),  where  X  is  the 
inverse  mean  free  path. 

Now,  in  the  theory  of  classical  liquids,  the  X-ray 
(neutron)  structure  factor  is  related  to  the  radial  distrib¬ 
ution  function  through  the  equation 

S(k)  =  1  +  p  |  [g(r)  -  Q  exp(i  k.r)  dr  (2.10) 

Various  theories  of  classical  liquids  relate  the  direct 
correlation  function,  C(r),  defined  by  the  equation 

C(r)  =  g(r)  (1  -  exp($(r)/kgT))  (2.11) 

to  the  effective  pair  potential,  $(r),  asymptotically 

C(r)  -  $ (r)/kgT  (2.12) 

By  the  method  of  molecular  dynamics  (Verlet  1967)  the 
pair  potential  obtained  from  the  liquid  structure  factor 
may  be  refined  so  that  it  reproduces  the  input  data.  Thus, 
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we  ,ire  able  to  determine  an  approximate  $(r)  from  the 
experimentally  observed  S(k). 

In  the  next  section  we  derive  expressions  for  the 
elastic  constants  and  phonon  frequencies  in  cubic  metals 
in  which  the  total  energy  is  given  by  equation  (2.8). 

The  short  range  repulsive  potential  is  examined  in  52.3 
and  values  for  the  repulsivity,  n,  are  estimated  for  the 
transition  metals.  Ni  and  Pt  are  examined  in  some  detail 
in  52.4  and  §2.5  and  central  pair  potentials  which  reproduce 
experimental  data  are  constructed.  Finally  in  52.6  we 
discuss  the  validity  of  central  pair  forces  in  the  f.c.c. 
metals  and  the  necessity  of  the  inclusion  of  Ey(V)  in  the 
expression  for  the  total  energy. 


52.2  Expressions  for  Elastic  Constants  and  Phonon 
Dispersion  Curves 

We  first  derive  expressions  for  the  elastic  constants 
of  a  cubic  crystal,  assuming  that  the  energy  is  given  by 
equation  (2.8),  by  the  method  of  homogeneous  deformation  of 
Fuchs  (1936).  We  then  derive  an  expression  for  the  dynam¬ 
ical  matrix  and  hence,  the  equations  governing  the  phonon 
dispersion  curves. 

a)  The  elastic  constants  of  a  cubic  crystal 

From  equation  (2.8)  the  energy  per  unit  undeformed 
volume  of  a  monatomic  crystal  can  be  written  as 

U  =  l  $(RS)  +  ev(V)  (2.13) 

o  s 

where  (Rs)  is  the  value  of  the  central  pair  potential  in 
equation  (2.9)  at  the  s  neighbour  shell.  is  the 

undeformed  atomic  volume  and 


Ey  (V)  =  Ey(V) 

"NV 

N  being  the  total  number  of  atoms. 


(2.14) 


The  elastic  constants  of  a  crystal  are  defined  as  the 


coefficients  in  the  expansion  in  powers  of  the  Lagrangian 
strain,  n^jj  of  the  change  in  the  internal  energy  density, 
due  to  a  homogeneous  deformation. 

Then  writing  the  change  in  the  internal  energy 
density  in  the  form 


AU  =  nijCij  +  *  nijnklCijkl 


(2.15) 


+  n  •  -  n,  ,ri  C..,,  + 

o  ij  kl  mn  ljklmn 


the  coefficients  C^j  ,  and  ^ijklmnare 

second  and  third  order  elastic  constants  respectively,  and 

are  evaluated  at  zero  strain. 


But,  in  the  approximation  of  equation  (2.8),  the 
change  in  internal  energy  can  be  written  as 


■(' 


4U  ■  'Bij  +  e’{i5  Hi 


*  ^  Bijkl  +  (ev  +  ev)6ij6kl  "  2Ev6jk6il)  n--ni 


Jijkl  •  (lv  '  tvJuijukl  ‘"vujkuilj  "ij "kl 

7  (Bijklmn  +  (*i  +  3ev  +Ev  j 6ij 6kl6r 
-  6  f  e_*.  +  e”  j  6 


mn 


il  jk  mn 


*  8  ci  4in8 jk4lm]  "ij"klnmn 


where 


R  -  *1  v  1  dji 

Bij  ~  2r  i  r  dr 

o  s 


r|  rJ 


B 


ijkl 


=  1  Y  1  d  f 1  dt] I  Rs  RS  RS  Rs 

2fl  *•  r  dr  Ir  dr  I  I  l  3  k  1 

os 


Bijklmn 


_  1  rid  fl  d  (1  d^ll 

‘  r  Hr  [r  Hr(r  drJJ 


(2.16) 


S( 2 . 17 


RS  pS  j^S  RS  j^8 


1  j  k  1  m  n 


and 


e(n)  _  yn  dnev 


dV 


n 


(2.18) 
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Comparing  equations  (2.15)  and  (2.16)  and  taking 
account  of  the  symmetry  conditions,  we  find  that  the 
second  order  elastic  constants  are  given  by  the  equations: 

CH  =  4-  I  (rs)(xs)4  +  e»  -  C'  (2.19) 

o  s 

c12  =  2^“  I  ♦"  (rs)(xs)2(ys)2  +  e"  +  (2.20) 

os 

C44  =  2S“  I  ♦  "  (rs)(xs)2(ys)2  -  (2.21) 

o  s 

with  the  equilibrium  condition 

C1  =  25“  l  *’<rs)(xS)2  +  e '  (2.22) 

o  s 

where 

. (n),  s>  fl  d  In  ^ ,  s. 

'■rs  drs-> 

The  Cauchy  discrepancy  Martin  (1975)  is  then 

C12  -  C44  *  *  2cv  (2-23> 

We  note  that  for  ev  =  0,  i.e.,  for  central  forces, 
this  discrepancy  is  zero,  which  is  just  what  we  would  expect. 

The  equations  for  the  third  order  elastic  constants 
follow  naturally 

Clll  ■  25—  l  *"'<rB>-<xB>6  4  3c;  -  3c;  .  e;" 

o  s 

C112  =  75"  l  +’"(rs)-(xs),,<ys)2  -c;  4  c;-4c"; 
o  s 

C166  *  75-  l  *"'<rS)-(xE)4(yS)2  4  c;  -  c; 
o  s 

c144  *  75-  l  *'"<rB)>(xs)2(ys)2(zs)2  -  c;  -  3c; 
o  s 

C123  =  2 ir  l  ^M,(rs)*(xs)2(ys)2(zs)2  +  +  3e  "  +  e»’ 

o  s 

C4S6  *  75-  l  *'"<rB>-<xB)2<yB>2<zB>2  4 
o  s 


A  little  algebraic  manipulation  leads  to 


'12 


-  C44  =  C 


456 


-  C 


144 


=  C 


166 


"  C112  +  C123 


-  C 


456 


=  2e^  +  e ”  (2.25) 

which  is  a  generalised  Cauchy  condition  for  the  model 
described  by  equation  (2.8).  This  condition  provides  a 
measure  of  the  suitability  of  the  approximation  made  in 
equation  (2.8)  for  particular  cubic  metals. 


b)  Phonon  Frequencies 

In  the  harmonic,  and  pair  potential  approximation  the 
phonon  frequencies,  <u(k)  are  given  by  the  equation 


u>  (k)  U  (k)  =  T  D  0  (tcx’lk)  U  (k ' )  (2.26) 

O  KTg  “P  —  P 

where  the  dynamical  matrix  D(k)  is  defined  by  the  relation 
Dag(<K'|k)  =  l  Wap  exp  (-ik*(x(o)  -x(l’))) 


[W] 


i 


(2.27) 


with 


(2.28) 


W  is  the  potential  energy  of  the  crystal  and  Ua  is  the 

displacement  is  the  a  direction  of  the  ic**1  basis  atom,  of 

mass  M  ,  in  the  l**1  unit  cell. 
k  ’ 


From  equations  (2.27)  and  (2.29)  it  follows  that 

^  .i  (si  b  •  •- 


- 11  -  *01-  •«')  ♦.*  1°  i;) 


x  exp  (-  i  k  «(x(o)  -  x  (I*)))  + 


♦ 


«  «(e  1  ♦ 
oB  v 


where  the  asterisk  denotes  the  sum  is  only  taken  o vet* 
and  f°|  being  nearest  neighbours. 


§2.3  The  short  range  potential 


We  showed  in  the  introduction  to  this  chapter  that 
the  short  range  repulsive  potential  between  two  ions  when 
they  are  brought  close  together  can  be  written  as 


B_ 

♦  (r)  n  =  A  exp  ( -  br) 

r  r  e 


(2.30) 


By  studying  various  experimental  data,  we  estimate  the 
value  of  n  for  some  f.c.c.  transition  metals. 

a)  Radiation  Damage 

Various  workers  have  studied  the  displacement  threshold 
energies  in  radiation  damage  and  the  work  has  been  exten¬ 
sively  reviewed  by  Vajda  (1977).  The  measured  threshold 
energies  of  a  number  of  metals  have  been  found  to  be  consis¬ 
tent  with  the  parameter  b  in  equation  (2.6),  taking  values 
between  4  and  6  A  .  Anderson  and  Sigmund  (1965) 

approximately  represented  the  short  range  interaction  bet¬ 
ween  two  atoms  in  the  metal  by  the  form 


,(r)  =  A0Z^^  exp  (-br) 


(2.31) 


where  Z  is  the  atomic  number  and  AQ  and  b  have  values  52eV 
and  4.56  A°-*  respectively. 


The  melting  temperature  at  zero  external  pressure 


The  simplest  and  most  successful  theory  of  melting  is 
that  of  Lindemann  (1910).  In  this  theory,  the  solid-liquid 
transition  occurs  when  the  root  mean  square  displacement  of 


each  lattice  atom  due  to  thermal  agitation,  becomes 
about  a  tenth  of  the  nearest  neighbour  distance. 


All  theories  of  melting  which  take  into  account  the 
potential  between  atoms  show  that  the  melting  temperature 
Tmeit»  is  proportional  to  the  strength  of  the  potential 
(Hoover  and  Ross,  1971).  Thus,  for  the  inverse  power 
potential  given  in  equation  (2.30) 


T 


melt 


where  a  is  some  constant. 


(2.32) 


In  order  to  relate  this  to  the  Born-Mayer  potential  we 
first  equate  the  two  forms  of  the  potential  in  equation 
(2.30)  in  the  energy  range  over  which  they  are  valid,  we 
find  that  for  the  f.c.c.  metals. 


and 


n  =  (_/S-l)  ab 
in3 

e  =  A  exp  (-ab/2) 


(2.33) 


(2.34) 


Thus,  the  melting  temperature  is  related  to  the  Born- 
May  parameter,  b,  by  the  equation 


T 

melt  _  y*  exp  (-ab/2) 

~wr  -  i - 

(0.9/2)(/I'1>  ab/ln3 


(2.35) 


where  y ^  is  some  constant. 

3/2 

In  figure  2.1  we  have  plotted  Tmelt/Z  against  ab 
from  equation  (2.35)  for  some  f.c.c.  transition  metals, 
with  b  taken  to  be  4.56  A°”*.  The  small  deviation  of  the 
experimental  points  from  the  curve  indicates  the  usefulness 
of  the  forms  of  the  potential  given  in  equation  (2.30) 


Fuchs  (1936)  and  Huntington  (1953),  who  used  the  Born- 
Mayer  form  of  the  potential  to  calculate  the  elastic  constants 
of  Cu,  estimated  the  parameter  b  to  be  between  6  and  7. 

This  results  in  n  taking  values  between  14  and  17.  Radiation 
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damage  analysis  of  the  transition  metals  (Vajda  1977) 
also  leads  to  values  of  b  consistent  with  the  above.  The 
hardness  of  the  core  of  the  transition  metal  can  also  be 
found  if  we  examine  the  data  of  melting  temperatures  at 
high  pressures. 

c)  Melting  curves  under  pressure 

It  can  be  shown  (Hoover  and  Ross  1971)  that  for  purely 
repulsive  potentials  of  the  form  Vr11  the  melting  curve  is 
consistent  with  a  Simon  coefficient  (Simon  1937),  c,  given 
by 

c  =  1  +  3/n  (2.36) 

The  melting  data  of  the  transition  metals  (see  for 
example  Cannon  1974)  shows  that  c  is  only  very  slightly 
larger  than  unity,  implying  a  large  value  of  n. 

Thus  all  the  experimental  evidence  shows  that  n,  for 
the  transition  metals  is  of  the  order  of  16. 


12.4  Pair  potentials  and  Phonon  Frequencies  in  Ni 

Johnson  et  al.  (1976)  obtained  a  pair  potential  for 
liquid  Ni  from  the  measured  neutron  scattering  structure 
factor  data.  The  potential,  they  reported,  was  modified 
slightly  so  that  it  reproduced  the  experimentally  determined 
radial  distribution  function  in  a  molecular  dynamics  cal¬ 
culation.  We  have  used  this  as  a  starting  point  in  our 
attempt  to  calculate  a  central  pair  potential  for  Ni. 

It  is,  of  course,  recognised  that  in  Ni  the  force 
fields  may  well  be  different  in  the  liquid  and  solid  states 
since : 

a)  There  is  a  volume  change  on  melting  and  one  expects 
from  the  basic  theory  that  interionic  potentials  will  be 
density  dependent.  This  is  a  direct  consequence  of  the 
definition  of  the  pair  potential. 

b)  The  Fermi  surface,  enclosing  the  electrons  which 
screen  the  Ni  ions  in  the  metal ,  has  a  complicated  shape  in 
the  crystal.  One  would  expect,  in  contrast,  the  Fermi 
surface  in  the  liquid  to  be  isotropic.  Naturally,  on  melting 


detailed  changes  will,  therefore,  occur  in  the  screening 
and  hence,  in  the  effective  interionic  potential. 

c)  The  electronic  mean  free  path  is  shorter  in  the 
liquid  than  in  the  solid.  Oscillations,  due  to  the  incom¬ 
plete  screening  of  ions,  will  therefore,  from  the  Uncer¬ 
tainty  Principle,  be  more  heavily  damped  in  the  liquid 
than  in  the  solid. 

In  view  of  our  present  lack  of  knowledge  of  the  inter¬ 
ionic  forces  in  the  transition  metals,  we  have  assumed  that 
the  changes  due  to  (a)  and  (b)  above,  are  hardly  likely  to 
be  larger  than  the  uncertainties  in  such  potentials.  The 
effect  of  the  greater  mean  free  path  in  the  liquid  can  be 
accounted  for  by  multiplying  the  effective  interionic 
potential  by  an  exponential  damping  factor. 

We  have  studied,  therefore,  how  well 
(i)  the  elastic  constants 
(ii)  the  phonons  throughout  the  B7. 
can  be  reproduced,  by  making  small  changes  in  the  potential 
extracted  from  the  liquid. 

From  the  arguments  given  in  §2.1  we  expect  the  central 
pair  potential  in  Ni  to  have  the  form  given  by  equation 
(2.9) 

$(r)  =  ^  cos  (2kp,r  +  e) 

r  r 

To  take  account  of  the  electronic  mean  free  path  this 
equation  can  be  modified  to  read 

$(r)  =  ^  cos  (2k„r  +  6)  exp  (-Xr)  (2.37) 

rn  r  * 

where  X  is  the  inverse  of  mean  free  path 

From  §2.3  we  argued  that  n  was  of  the  order  of  16. 

To  obtain  the  other  parameters  in  the  potential  we  have 
relied  on  the  liquid  pair  potential  of  Johnson  et  al.  (1976). 
The  values  of  kp  and  6  were  estimated  from  the  wavelength 
and  nodal  positions  of  the  long  range  part  of  the  liquid 
potential.  Neglecting  X,  A  and  B  were  determined  by  fitting 


the  position  and  depth  of  the  first  minimum  of  this  poten¬ 
tial.  These  values  are  tabulated  in  table  2.1. 

The  mean  free  path,  L,  in  liquids  is  about  twice  the 
nearest  neighbour  distance.  In  solids  the  mean  free  path 
is  much  greater,  and  we  have  assumed  the  ratio  of  the  mean 
free  path  in  the  liquid  to  that  in  the  solid,  to  be  of  the 
order  of  the  ratio  of  the  room  temperature  to  the  melting 
temperature.  Hence,  Lson<j  *  5  This  gives  X  a 

value  of  about  0.0 5A°-^  in  solid  Ni. 

The  contribution  of  this  potential  to  the  phonon  freq¬ 
uencies  was  calculated  by  the  use  of  equations  (2.26)  to 
(2.29).  The  potential  was  truncated  at  the  fourth  neigh¬ 
bour  distance  and  X  was  taken  to  be  zero.  We  have  in  fact 
calculated  w(k)  at  different  cut-offs  (up  to  35.2A°)  with 
a  finite  value  of  X  and  the  discrepancies  between  the  res¬ 
ults  were  found  to  be  less  than  5%,  well  within  the  range 
of  experimental  uncertainties.  However,  we  find  that  with 
this  potential,  X(k)  in  some  modes  turn  out  to  be  imag¬ 
inary,  indicating  some  kind  of  instability  in  the  solid 
when  phonons  are  propagated  along  these  directions.  To 
remove  ties e  instabilities,  we  have. varied  A,  B  and  6  with 
the  constraint  that  the  depth  of  the  first  minimum  is  kept 
roughly  the  same.  The  new  values  of  the  parameters  are 
also  tabulated  in  table  2.1  under  the  column  entitled 
'solid'.  The  depth  of  the  first  minimum  was  kept  the  same 
because  it  is  related  to  the  height  of  the  first  peak  in 
the  g(r)  of  the  liquid,  which  we  require  to  be  the  same. 
Both  these  potentials  are  drawn  in  figure,  2. 2.  It  is  inter 
esting  to  note  that  the  'solid  potential'  has  its  minimum 
displaced  slightly.  A  simple  way  in  which  the  potential 
can  be  made  density  dependent  is  by  scaling  the  interatomic 
distances.  The  volume  change  on  melting  would,  thus, 
result  in  the  distance  of  the  first  minimum  in  the  liquid 
potential  being  fractionally  greater  than  that  of  the 
'solid  potential',  which  is  consistent  with  our  potentials 
for  Ni, 


TABLE  2.1 


Parameters  of  the  solid  and  liquid  potentials  in  Ni. 


PARAMETER 

LIQUID 

SOLID 

A  (eV  X  16) 

44351.0 

40319.0 

B  (eV  X  3) 

2.26 

2.06 

2kf  (X-1) 

3.25 

3.25 

e 

0.994 

1.600 
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Having  obtained  a  reasonable  central  pair  potential, 
phonon  dispersion  curves  were  recalculated  with  the  inclu¬ 
sion  of  the  volume  dependent  energy.  Ev(V)  was  represented 
by  the  semi-empirical  form 


E  (V)  =  d 
V  ^71 


(2.38) 


which  was  prompted  by  the  free  electron  volume  dependence. 

Of  course,  this  does  not  have  a  first-principle  justific¬ 
ation  for  the  transition  metals,  but  it  appears  adequate 
for  our  present  purposes.  The  parameter  d  was  determined 
from  the  equilibrium  condition  expressed  by  equation  (2.22). 

We  have  plotted  in  figure  2.3  for  Ni  the  calculated 
w(k),  with  and  without  the  inclusion  of  Ev,  versus  k  in 
the  three  symmetry  directions  together  with  the  experimental 
points  (Birgeneau  et  al.  1964).  The  dashed  lines  are  those 
curves  calculated  from  the  central  pair  potential  of  the 
solid  and  the  full  lines  are  those  obtained  with  the 
inclusion  of  Ey.  It  may  be  observed  that  the  term,  Ev, 
provides  a  necessary  contribution  to  u)(k),  to  bring  it  in 
agreement  with  the  experimental  data. 

From  the  nature  of  E(V)  discussed  in  52.1,  we  see  that 
d  in  equation  (2.38)  must  be  positive.  The  actual  form  may 
not  be  unique,  but  the  representation  of  Ev(V)  by  a  single 
volume  dependent  term  is  the  most  useful,  as  the  one  para¬ 
meter  is  fixed  by  the  equilibrium  condition.  We  have 
repeated  the  calculations  for  different  forms  of  Ev, 
namely, 


Ev(V)  =  d/Vm 


where  m  was  of  the  order  of  unity. 


(2.39) 


The  changes  in  the  phonon  frequencies  from  the  m  =  2/3 
case  proved  to  be  quite  insignificant. 

We  have  also  calculated  the  second  order  elastic 
constants  of  Ni  with  and  without  the  inclusion  of  Ev(V), 
and  these  are  shown  in  table  2.2,  together  with  the  experi- 


22 


TABLE  2.2 


Second  order  elastic  constants  and 
bulk  modulus  B  for  f.c.c.  Ni. 


<1010N/m2) 

E(V)  -  0 

E(V)  =  d 

r m 

Expt. 

C11 

18.9 

24.1 

24.6 

C12 

11.7 

14.1 

15.0 

C44 

11.7 

13.2 

12.2 

B 

14.1 

17.4 

18.6 

(a)  de  Klerk  (1959) 
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mental  values  (de  Klerk  1958).  The  positive  d  in  equation 
(2.38)  gives  a  positive  Cauchy  discrepancy  in  accordance 
with  experiment. 

The  radial  distribution  function  for  Ni  using  the 
'solid  potential'  was  generated  by  a  molecular  dynamics 
calculation.  The  resulting  g(r)  does  not  differ  signifi¬ 
cantly  from  that  obtained  from  the  liquid  potential. 


S2.5  Pair  Potentials,  Phonon  Frequencies  and  Liquid 
Structure  Factor  in  Pt 

Recent  neutron  inelastic  scattering  measurements  from 
H  adsorbed  on  a  Pt  surface  by  Wright  et  al.  (1976)  have 
been  interpreted  in  general  terms  by  Mahanty  et  al.  (1976). 

It  is  clear  from  their  work  that  a  quantitative  interpret¬ 
ation  of  the  neutron  results  will  require  knowledge  of  the 
long  range  interionic  potential  for  Pt.  We  have,  therefore, 
attempted  to  construct  an  interionic  potential  for  Pt. 

The  arguments  presented  in  52.3  which  led  us  to  choose 
the  repulsivity  of  Ni  to  be  16,  also  apply  to  Pt.  This  is 
especially  true  of  the  data  on  the  melting  curves  under 
pressure.  With  this  value  of  n,  the  parameters  A  and  B  were 
adjusted,  so  that  the  first  minimum  of  the  potential  appeared 
at  about  the  nearest  neighbour  position  in  Pt  with  depth 
approximately  equal  to  that  in  Ni.  An  initial  value  for  2kp 
was  estimated  from  the  value  for  Ni,  but  with  the  increased 
atomic  volume  taken  into  account.  The  phonon  frequencies 
in  the  three  symmetry  directions  were  calculated  with  this 
potential.  The  values  of  2kp  and  e  were  then  varied  until 
a  good  fit  to  the  phonon  dispersion  curves  was  obtained. 

We  were,  thus,  able  to  find  a  useful  representation  of  the 
central  pair  potential  in  Pt,  namely. 

0 (r )  =  409732  +  3.29  cos  (3.49r  +  1.55)  (2.40) 


where  r  is  measured  in  A°,  and  the  potential  is  truncated 
at  the  fourth  neighbour  distance.  This  potential  is 
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plotted  in  figure  2.4. 

Then,  including  Ev(V)  in  the  form  given  by  equation 
(2.38),  the  phonon  frequencies  were  recalculated  and  the 
results  are  shown  in  figure  2.5,  together  with  the  experi¬ 
mentally  measured  values  (Dutton  et  al.  1972).  Unlike  in 
the  case  of  Ni,  the  effect  of  Ey  turns  out  to  be  quite 
small. 

The  general  agreement  between  the  calculated  and 
experimental  values,  with  the  exception  of  the  pronounced 
Kohn  anomaly  in  the  (110)  direction,  is  very  good.  It  may 
be  that  to  reproduce  the  Kohn  anomaly,  one  would  need  to 
take  account  of  the  long  range  oscillatory  tail.  However, 
if  this  were  the  case,  Kohn  anomalies  would  also  appear  in 
the  theoretical  calculations  of  the  dispersion  curves  of  Ni. 

That  this  is  not  experimentally  observed  leads  one  to  think 
that  it  is  the  detail  of  the  Fermi  surface  that  is  important 
(see  for  example  Flores  et  al.  1979). 

Recently,  Tripathi  and  Nand  (1979)  calculated  the  phonon 
frequencies  of  Pt  using  a  model  potential  approach.  Though 
their  results  are  in  fair  agreement  with  experiment,  they  were  unabb 
to  reproduce  the  Kohn  anomaly.  This  is  in  spite  of  the  fact 
thal  their  calculations  were  made  in  reciprocal  space  and 
so  took  account  of  the  long  range  tail.  Their  estimates  of 
the  magnitude  of  the  anomaly  turns  out  to  be  very  small. 

Their  calculations,  thus,  show  that  the  inclusion  of  the 
long-range  tail  does  not  necessarily  result  in  observable 
Kohn  anomalies  and  the  observation  of  such  should  be  attri¬ 
buted  to  the  details  of  the  Fermi  surface. 


Wc  have  also  calculated  the  second  order  elastic  cons¬ 
tants  of  Pt  and  these  are  tabulated  together  with  the  experi¬ 
mental  values  (Gschneider  1964)  in  table  2.3.  Again,  we  note 
that  Ey  does  not  appear  to  make  any  significant  contribution. 
The  truncation  of  the  interatomic  potential  is  not  a  good 
approximation  in  the  evaluation  of  the  elastic  constants, 
because  from  equations  (2.19)  to  (2.21)  we  see  that  the 
expressions  for  them  are  not  absolutely  convergent  if  the 
mean  free  path,  1/X,  is  infinite. 


TABLE  2.3 


Second  order  elastic  constants  and 
bulk  modulus  B  for  f.c.c.  Pt. 


(1010N/m2) 

E( V)  =  0 

E(v)  ■ 

Expt . 

C11 

22.8 

23.4 

C1 2 

9.9 

10.2 

C44 

9.9 

10.0 

B 

14.2 

15.6 

27.8 

(a)  Gschneider  (1964) 
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In  contrast  with  our  work  on  Ni,  we  decided  to  test 
our  potential  in  Pt,  obtained  by  fitting  the  phonon  freq¬ 
uencies,  by  calculating  the  liquid  S(k)  and  comparing  it 
with  experiment. 

From  equation  (2.40)  we  see  that  the  core  is  very 
hard,  so  we  expect  S(k)  to  have  the  main  features  of  the 
known  hard  sphere  form.  This  latter  result  is  independent 
of  temperature,  depending  only  on  the  ratio  of  the  volume 
occupied  by  the  hard  spheres  to  the  total  volume  (Ashcroft 
and  Lekner  1966).  We,  therefore,  ran  a  moleculer  dynamics 
program  at  a  higher  temperature  than  the  known  melting 
temperature  of  Pt ,  to  ensure  that  the  solid  had  melted. 

In  the  program,  the  system  consisted  of  256  particles 

interacting  between  themselves  through  the  central  pair 

potential  given  by  equation  (2.40).  The  system  was  allowed 

—15 

to  evolve  for  five  hundred  time  steps,  each  of  10  seconds. 
The  time  interval  was  chosen  such  that  the  total  energy  of 
the  system  was  conserved  to  four  significant  figures.  The 
mean  temperature  was  then  calculated  over  the  next  two 
hundred  time  step  intervals  to  verify  that  equilibrium  had 
indeed  b^en  reached.  At  the  end  of  each  time  interval  the 
g(r)  and  its  fourier  transform,  S(k)  were  calculated  accor¬ 
ding  to  equation  (2.10).  There  was  no  significant  differ¬ 
ence  at  each  time  interval. 

We  have  plotted  S(k)  versus  k  for  Pt  in  figure  2.6. 

The  crosses  indicate  the  positions  of  the  maxima  and  minima 
of  the  experimental  results  of  Waseda  and  Ohtani  (1975a). 

The  close  agreement  between  the  results  indicates  that 
the  potential  given  by  equation  (2.40)  also  provides  a  good 
description  of  the  liquid  state. 

Another  interesting  observation  is  that  the  agreement 
is  very  good  even  at  the  third  peak.  The  long  range  part 
of  S(k)  is  dependent  on  the  repulsivity,  therefore,  this 
agreement  is  a  vindication  of  our  choice  of  n  for  Pt. 


12.6  Discussion 


In  the  above  sections  we  have  obtained  central  pair 
potentials  for  Ni  and  Pt .  The  usefulness  of  these  pair 
potentials  depends  on  the  validity  of  the  approximation 
expressed  by  equation  (2.8).  We  have  also  assumed  that 
the  central  pair  forces  in  the  solid  and  liquid  states 
are  very  nearly  the  same.  We  now  consider  whether  these 
approximations  are  good  in  the  f.c.c.  transition  metals. 

In  12.2  we  saw  that  if  the  approximation  represented 
by  equation  (2.8)  is  to  be  valid,  in  any  cubic  metal,  then 
their  second  and  third  order  elastic  constants  must  satisfy 
equation  (2.25).  Since  this  work  was  first  completed 
(Matthai  1978a)  Cousins  and  Martin  (1978)  showed,  by  con¬ 
sidering  inner  displacements  in  a  crystal,  that  the  linear 
combination  of  elastic  constants 

Q3  =  C12  "  C44  +  C144  "  C456  (2.41) 

and 

Q6  =  2(C12  "  C44}  _  C123  +  C144  +  C112  "  C166  (2*42) 

are  devoid  of  contributions  from  central  or  isotropic  inter¬ 
actions.  This  is  completely  consistent  with  our  equation 
(2.25) 

We  have  tabulated  in  table  2.4,  the  three  linear  com¬ 
binations  of  elastic  constants  appearing  in  equation  (2.25) 
as  dete..  ..lined  from  experiment  for  some  f.c.c.  metals  for 
which  the  third  order  elastic  constants  have  been  measured 
and  for  directionally  bonded  diamond-type  Si  and  Ge. 

Although  these  conditions  are  not  satisfied  by  any  of  the 
solids  considered,  it  appears  that  it  is  a  better  approxim¬ 
ation  in  the  f.c.c.  metals  than  in  Si  or  Ge.  It  should  be 
noted  that  the  errors  in  the  measurements  of  the  third 
order  elastic  constants  are  quite  large. 


TABLE 


(c)  Thomas  (1968) 

(d)  McSkimmin  and  Andreatch  (1964) 

(e)  Drabble  and  Fendly  (1967) 
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Vie  have  also  calculated  the  third  order  elastic  cons¬ 
tants  of  Ni  and  they  are  tabulated  with  the  experimental 
values  in  table  2.5.  With  the  exception  of  C^2  3»  the  agree¬ 
ment  between  the  two  sets  of  values  is  very  good,  especially 
if  one  recognises  that  to  obtain  the  third  order  elastic 
constants  one  requires  the  third  derivative  of  the  inter 
atomic  potential. 

This  gives  added  weight  to  our  assumptions  that  angu¬ 
larity  in  the  force  fields  is  less  prevalent  in  the  f.c.c. 
than  in  the  b.c.c.  transition  metals.  Thus,  it  would  seem 
that  with  the  phonon  frequencies  and  liquid  S(K)  well 
reproduced,  the  idea  of  central  pair  forces  in  the  f.c.c. 
transition  metals  is  a  useful  one. 

Another  aspect  of  our  work  has  been  the  explicit 
introduction  to  the  volume  dependent  term  in  the  total 
energy.  This  enabled  us  to  determine  its  contribution  to 
the  elastic  constants  and  the  phonon  frequencies.  It  also 
results  in  a  positive  Cauchy  discrepancy  for  isotropic  forces 
as  is  observed  in  most  metals.  Furthermore,  by  introducing 
the  parameter  d  in  equation  (2.38),  we  were  also  able  to 
satisfy  the  equilibrium  condition  given  by  equation  (2.22) 
From  our  calculations ,  it  turns  out  to  be  ir  jch  smaller  than 
the  coefficient  of  the  leading  term  in  the  free  electron  gas 
energy  expansion.  Previous  workers  were  not  able  to  satisfy 
the  correct  equilibrium  condition  and  get  reasonable  phonon 
frequencies  because  of  the  form  they  used  for  Ev.  A  review 
of  the  literature  shows  that  most  workers  have  used  incorrect 
equilibrium  conditions  with  E  as  the  free  electron  energy. 
Upadhyaya  (1978)  pointed  this  out,  but  suggested  that  to  get 
a  positive  Cauchy  discrepancy  one  needs  to  go  to  higher 
order  in  the  perturbation.  We  have  shown  that  this  need 
not  be  the  case  if  Ev  were  not  just  the  free  electron  energy 
and  that  part  arising  from  charge  neutrality,  but  also  con¬ 
sisted  of  other  many  body  terms,  which  are  averaged  out  to 
make  an  isotropic  contribution. 


TABLE 


Finally,  we  conclude  this  chapter  with  some  comments 
on  the  accuracy  of  the  pair  potentials  obtained.  Varia¬ 
tions  on  the  accuracy  of  the  parameters  A,  B,  k^  and  e 
result  in  changes  in  the  dispersion  curves.  In  so  far  as 
the  phonons  with  k  i  0  accurately  determine  a  potential, 
we  estimate  these  parameters  to  be  accurate  to  better  than 
5%.  However,  this  is  for  a  fixed  value  of  X.  Because  of 
the  truncation  of  our  potential,  a  further  uncertainty  due 
to  the  parameter  X  is  also  introduced.  So,  before  claims 
of  an  accurate  potential  can  be  made,  the  actual  value  of 
X  is  required. 
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CHAPTER  3 


The  angularity  of  the  charge  distribution 
(force  fields)  in  metals 

§3.1  Introduction 


From  linear  response  theory,  it  follows  that  the  force 
fields  in  metals  are  directly  related  to  the  charge  distri¬ 
bution.  Although  we  assumed  that  central  forces  provide  a 
good  description  of  the  f.c.c.  transition  metals,  and  showed 
that  the  assumption  to  be  a  good  one  in  calculating  the 
elastic  constants,  phonon  frequencies  and  the  liquid  struc¬ 
ture  factor;  we  realise  that  a  full  description  of  the 
force  fields  would  have  to  take  account  of  the  angularity 
in  the  charge  distribution.  This  angularity  is  most  easily 
discernible  in  the  b.c.c.  transition  metals.  The  inadequacy 
of  the  central  pair  potential  model  for  the  b.c.c.  transit¬ 
ion  metals  is  most  apparent  in  the  behaviour  of  their  elas¬ 
tic  constants  and  phonon  dispersion  curves. 

With  the  assumption  of  a  central  pair  potential  plus 
a  volume  dependent  energy,  we  saw  that  the  calculated  Cauchy 
discrepancy,  ~  ^44’  will  be  small  and  positive.  An 

examination  of  the  second  order  elastic  constants  of  the 
b.c.c.  transition  metals  reveals  that  their  Cauchy  discrep¬ 
ancies,  though  positive  (with  the  exception  of  Cr) ,  differ 
a  great  deal  in  magnitude.  Alone  among  the  metals,  Cr  dis¬ 
plays  a  negative  discrepancy,  a  feature  also  found  in  the 
covalently  bonded  solids,  Si  and  Ge.  From  our  formulation 
of  Ev,  this  can  never  arise  as  d  given  by  equation  (2.38) 
is  taken  to  be  positive  and  so  to  correct  this  deficiency 
in  the  model,  angular  forces  must  also  be  included  in  the 
expression  for  the  total  energy.  It  further  indicates  that 
perhaps  this  angularity  can  be  represented  in  terms  of  bond 
hybrids . 

Further  evidence  of  the  angularity  in  the  b.c.c.  trans¬ 
ition  metals,  from  an  observation  of  certain  linear  combin¬ 
ations  of  elastic  constants,  has  been  given  by  Chikazumi 


(1979).  By  appealing  directly  to  the  measured  second  order 
elastic  constants,  he  was  able  to  show  that  for  Ni  the 
atomic  force  is  almost  central,  wheras  he  concluded  that 
the  atomic  shape  of  the  Fe  atom  in  the  low  spin  state  is 
anisotropic . 

Brockhouse  et  al.  (1967)  made  a  Born- von  Karman  force 
constant  fit  to  the  phonon  dispersion  data  of  Fe,  and  found 
that  the  1XX  component  was  greater  than  the  1XY  component. 
They  interpreted  this  as  being  indicative  of  an  attractive 
force  between  nearest  neighbours.  This  observation  was  con¬ 
firmed  by  Minkiewicz  et  al.  (1967).  Brockhouse  et  al. also 
observed  that  the  1XX  component  was  greater  than  the  1XY 
component  for  all  the  b.c.c.  transition  metals  that  had 
been  investigated.  Interestingly,  they  found  the  opposite 
inequality,  which  implies  repulsive  forces  between  nearest 
neighbours,  as  would  be  intuitively  expected  to  be  true  in 
the  non-transition  b.c.c.  metals.  Analysis  of  all  the 
f.c.c.  metals  also  show  first  neighbour  repulsion.  These 
facts  rather  strongly  suggest  that  in  the  b.c.c.  transition 
metals,  the  d-electrons  are  involved  in  some  kind  of  quasi- 
covalent  bonding. 

Considerable  information  about  the  charge  density  can 
be  obtained  from  an  analysis  of  the  X-ray  scattering  factor 
data  of  solids.  A  useful  way  of  extracting  this  information 
is  to  compare  the  X-ray  scattering  results  from  the  solid 
state  with  the  free  atom  prediction.  This  tells  us  how  the 
basic  building  block  in  the  solid  suffers  from  that  in  the 
free  atomic  state. 

In  the  rare  gas  solids  (Ar,  Kr  and  Xe)  the  measured 
scattering  factors  show  no  difference  from  the  free  atom 
results,  which  is  in  agreement  with  the  theory  that  in  the 
solid  state  the  building  blocks  are  still  the  free  atoms. 

A  derivation  from  the  free  atom  results  is,  however,  obser¬ 
ved  in  the  simple  and  d-band  metals  and  the  covalently 
bonded  solids. 

In  covalently  bonded  Si,  once  angularity  is  put  into 
the  charge  distribution  around  the  ionic  centre  to  take 
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account  of  the  bonding,  the  discrepancy  between  the  scatt¬ 
ering  factors  in  negligible  (De  Marco  and  Weiss  1965). 

This  is  in  direct  contrast  with  the  metals,  where  the  dis¬ 
crepancy  persists  even  when  angularity  is  introduced  into 
the  ionic  centre  charge  distribution  (Prown  and  Taylor 
1969). 

Batterman  et  al.  (1960)  measured  the  X-ray  scattering 
factors  of  Fe,  Cu  and  Al.  Their  results  were  found  to  be 
about  4%  smaller  than  the  calculated  free  atom  scattering 
factors,  with  Fe  showing  the  greatest  discrepancy.  Brown 
and  Taylor  (1969)  suggested  that  this  discrepancy  is  a 
feature  of  metallic  bonding,  which  is  responsible  for  the 
modification  of  the  electronic  charge  distribution.  Exper¬ 
iments  by  Paakari  and  Suortti  (1967)  on  Fe ,  Hosoya  and 
Yamagishi  (1966)  on  Cu  and  Inkinen  and  Suortti  (1964)  on 
Ni  confirmed  the  discrepancy  between  the  measured  and  free 
atom  values  found  by  Batterman  et  al.,  although  they  found 
it  to  be  much  smaller.  In  Fe  the  discrepancy  was  found  to 
be  about  1.7%  and  extended  up  to  values  of  sin  e/X  =  0.5$  *. 
The  discrepancies  in  Cu  and  Ni  are  smaller  and  vanish  almost 
completely  at  the  second  reflection,  which  is  very  close  to 
the  first  one  in  angular  position.  «This  appears  to  indi¬ 
cate  that  the  electronic  redistribution  in  Fe  is  greater 
than  in  Cu  or  Ni. 

In  the  next  section  we  consider  the  theory  of  X-ray 
scattering  from  solids  using  a  Wannier  function  formalism. 

In  §3.3  we  use  the  results  obtained* from  this  approach  to 
interpret  the  experimental  scattering  data  of  Ihe  transit¬ 
ion  metals.  We  pay  special  attention  to  Fe  because  of  the 
availability  of  reliable  scattering  data  up  to  large  scatt¬ 
ering  angles.  We  discuss  the  nature  of  the  charge  distrib¬ 
ution  in  the  transition  metals  and  the  implications  on  the 
force  fields  in  §3.4 


§3.2  X-ray  Scattering  in  Solids  in  the  Wannier  Function 
Formalism 


The  X-ray  scattering  factor  S(k)  is  given  by 
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C_k>  =  |y  p(r)  exp  (ik.r)  dr 


(3.1) 


where  p (r )  is  the  charge  density  of  the  crystal  of  volume 
V. 


Let  us  consider  a  crystal  of  N  atoms  with  ion  core 
electronic  wavef unctions ,  ^c(r)  and  valence  electrons, 
*Pk(r).  The  total  charge  density  of  the  crystal  is 


p(r)  =  l  l  ip  (r-R  )  ip  (r-R  ) 
—  £  f*  c  —  — n  c  —  — n 


+  I  *k(r) 

occupied  —  — 

levels 


(3.2) 


which  we  can  write  as 


p(r)  =  £  o(r-R  )  +  £  tp.  ( r )  ip.(r)  (3.3) 

n  occupied  —  — 

levels 

where  R  is  a  lattice  vector. 

— n 

The  ^^(£)  are  the  Bloch  functions  (BF's)  generated  by 
the  one-bolly  periodic  potential  V(r).  This  potential  can 
be  written  in  the  form 

V  (r )  =  Vh(r)  +  {excM  (3.4) 

6p 

where  V  (r)  is  the  Hartree  potential  and  e  the  exchange 
h  x  c 

and  correlation  energy  which  is  a  functional  of  the  density. 

For  simplicity,  we  assume  the  valence  band  to  be  des¬ 
cribed  by  one  occupied  non-degenerate  band.  The  BF*s  can 
then  be  written  in  terms  of  Wannier  functions  (WF's), 
a (r-R  ) ,  name ly 


(3.6) 
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Substituting  equation  (3.5)  into  the  expression  for 
the  charge  density,  we  find 

p(r)  =  l  pk(r)  +  l  a(r-R)  (3.7) 

k<kF  —  n  n 

where 

p.(r)  =  N  1  7  exp(ik.(R  -R  ))  a(r-R  )  a(r-R  ) 
k  —  *•  —  — m  — n  - m - n 


and  kp  is  the  Fermi  wavevector  which  defines  the  boundary 
of  the  Fermi  surface. 


To  proceed  further,  we  assume  that  the  Fermi  surface 
is  spherical.  Then  it  can  be  shown  that  the  scattering 
factor  is  given  by 


S(k) 


i 


mn  kr| R  -R  | 

F 1  — m  — n 1 


a(r-R  )  a(r-R  )  exp(ik.r)dr 
—  — m  —  — n  —  —  — 


+ 


1  o(r-Rn)  exp(ik.r)  dr 
n 


(3.8) 


Then  writing  R  -R  =  R  and  letting 
— m  — n  — 


r  ■>  r  +  Rm  +  R_ 
—  —  — m  — n 


in  the  first  term  in  equation  (3.8)  we  find 

SOc)  =  y  exp  [ii-'Em  *  S„>]  T(k) 

-  k-R  l  - 5 -  J  B.  - 


kpR 


+  l  expdk.R^)  fc(k)  (3.9) 

n 

where 

IR  =  [  a(r-R)  a(r+R)  exp(ik.r)dr  (3.10) 

"  J°  2  7 


I] 


f  (k)  =  c(r)  exp  (ik.r)dr 
c  ~  -  -  - 


(3.11) 


and  fl  is  the  atomic  volume.  Since, 


l  exp(ik.Rn)  =  )>  exp(iG.R^)  6(k-G)  (3.12) 


where  G  are  the  reciprocal  lattice  vectors  and  R.  is  the 

"til 

position  of  the  i  basis  atom  in  a  unit  cell,  the  scatt¬ 
ering  factors  are  only  defined  at  the  reciprocal  lattice 
vectors : 


KG)  =  Sq(G)  J  exp  -  -i 


1-’— nl  3^l(kFRn)  I„  (G)  +  f  (G) 

2  kvR  n 

r  n  — 1 


(3.13) 

where  the  geometrical  structure  factors  of  the  lattice  are 
defined  by 


5q(G)=  £  exp(iG.R^) 


(3.14) 


i\ow  we  can  write 


I  =  I 


n  R  i( equal  R  ) 
n  n 


where  . ,  %  _  N  is  that  part  of  £  such  that  |R  I  is  a 

i( equal  RnX.  *  n  1  n1 

constant  and  l  is  the  sum  over  all  shells  of  atoms  dis- 

Rn 

tance  Rn  from  some  central  atom.  Then,  defining  I^n(G)  to 
be  the  average  of  Ir^(G)  over  all  R^  with  |R^|  =  Rn,  we  can 
write  equation  (3.13)  in  the  form 


such  that  R  is  a 
1  n 1 


S(G)  =  [fc(G)  +  f0(G T\  5cCG) 


l  f„(G)  S„<G)  (3.15) 

all  superlattices 
labelled  by  n 


where 


S  (G)  =  l  exp  f1-'  f^i  +  -n] 
n"  ni  ll  2-JJ 


(3.16) 


define  the  geometrical  structure  factors  of  the  n  super- 
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lattice,  and 


f  (G) 
n  — 


kFRn 


ht  <S> 


(3.17) 


are  the  form  factors  of  the  scattering  centres  of  the  n^"*1 
superlatt ice . 


In  particular 


f  (G)  = 
o  — 


j)a(r)]^  exp(iG.r)dr 


(3.18) 


and  this  together  with  fQ(G)  make  up  the  form  factor  of  a 
single  ion.  The  position  vectors  of  the  scattering  centres 
of  the  superlattices  are,  from  equation  (3.16),  just  the  mid 
points  of  the  lines  joining  an  atom  with  its  neighbours. 


Hence,  once  the  WF’s  are  known,  the  X-ray  scattering 
factors  can  be  determined  by  use  of  equation  (3.15).  The 
importance  of  this  equation  is  that  it  shows  us  that  one 
centre  charge  distributions  alone  do  not  tell  the  whole 
story.  For  a  full  picture  the  two  centre  terms  expressed 
by  equation  (3.10)  have  also  to  be  considered. 


S3. 3  Interpretation  of  the  X-ray  Scattering  Data  in  Solids 

In  the  rare  gas  solids  there  are  no  valence  electrons 
and  hence,  all  the  scattering  comes  from  the  core  electrons. 
Also,  since  there  is  little  or  no  overlap  between  the  core 
electrons,  the  WF's  in  the  solid  will  be  very  nearly  the 
atomic  functions.  Hence,  the  calculated  free  atom  scatter¬ 
ing  factors  should  agree  with  the  experimental  values. 

Insulators  and  semi-conductors  are  characterised  by  a 
lack  of  conduction  electrons.  This  corresponds  to  a  filling 
up  of  the  whole  BZ.  We,  therefore,  need  not  make  the  approx 
imation  of  a  spherical  Fermi  surface  as  we  can  sum  over  the 
occupied  states  exactly,  namely 

T  exp(ifc.(R  -R  ))  =  I  exp(ik.(R  -R  ))  =  N  6(R  -R  ) 
lc«kj.  Ik  over  BZ 


(3.19) 


-  M2  - 


Hence, all  the  superlattice  terms  in  equation  (3.15)  will 
vanish  leaving  only  the  form  factor  of  the  ionic  charge 
distribution.  However,  this  one  centre  form  factor  will  be 
in  terms  of  the  WF's,  which  because  of  overlap  differ  quite 
considerably  from  the  free  atomic  wavef unctions . 

In  the  free  atomic  state  all  degenerate  levels  are 
filled  up  equally  and  by  Unsolds  theorem  the  resulting 
charge  density  is  spherically  symmetric.  The  orthogonality 
of  the  WF's  ensures  that  they  have  the  symmetry  of  the  latt¬ 
ice  and  ,  thus,  angularity  in  the  charge  distribution,  con¬ 
sistent  with  the  crystal  structure,  is  contained  in  the  WF's. 
This  explains  why  De  Marco  and  Weiss  (1965),  by  including 
angularity  in  the  one  centre  charge  distributions  to  account 

O 

for  sp  hybridization,  were  able  to  obtain  agreement  with 
the  experimentally  measured  scattering  factors  of  Si. 

In  the  metals,  because  the  valence  electrons  spill 
outside  the  BZ ,  the  two  centre  terms  in  equation  (3.15) 
are  important.  Hence,  any  calculation  which  considers  only 
one-centre  terms  will  not  be  in  agreement  with  experiment. 
Though  the  sum 


l  exp(i_k.  (Rm-Rn)) 

k<  kp 

cannot  be  calculated  exactly,  as  shown  above,  the  scatter¬ 
ing  factors  can  be  approximately  determined  if 

(i)  the  Fermi  surface  can  be  approximated  by  a 
spherical  surface. 

(ii)  the  WF's  are  known. 

We  have  made  such  a  calculation  for  Be  and  this  is  presented 
in  the  next  chapter. 

The  difficulty  in  obtaining  the  WF's  in  metals  promp¬ 
ted  us  to  simplify  equation  (3.15)  into  a  useful  form,  so 
that  we  may  obtain  information  on  the  charge  density  in  the 
transition  metals  by  examining  the  experimentally  determined 
X-ray  scattering  factors. 


Stenhouse  et  al.  (1977)  calculated  the  scattering  fac 
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tors  of  Si  using  a  bond  charge  model.  In  this  model,  spher¬ 
ical  charge  distributions  were  placed  at  the  ionic  centres 
and  midway  between  nearest  neighbours  representing  the 
bonding  charge,  such  that  the  charge  density  in  LCAO  theory 
was  well  reproduced.  The  success  of  this  model  and  that  of 
De  Marco  and  Weiss  (1965)  suggests  that  the  angularity  in 
the  one  centre  terms  can  be  represented  by  bond  charges. 
Hence,  the  first  term  in  equation  (3.15)  can  be  written  in 
the  form 

f  (G)  5(G)  =  f  °(G)  S  (G) 
o—  o—  c  —  o  — 

+  l  f  b(G)  Sn(^)  (3.20) 

all  superlattices  n 
labelled  by  n 

c  b 

where  f  (G)  and  f  (G)  are  the  form  factors  of  the  charges 
(which  need  not  be  spherical)  at  the  ion  and  bond  sites 
respectively. 

Then  substituting  equation  (3.20)  into  equation  (3.15) 
we  obtain  an  approximate  formula  for  the  scattering  factor: 

S(G)  =  F  (G)  5(G)  +  7  F  (G)s  (G)  (3.21) 

0  o  L ,  ...  n  —  n  — 

all  superlattxces 

labelled  by  n 

where 

F0(G)  =  fc(G)  +  fQC(G) 

is  the  form  factor  of  the  ionic  distribution,  and 

F  (G )  =  f b(G)  +  f  (G) 
n  —  n  —  n  — 

represent  the  form  factors  of  the  bond  charges. 

The  size  of  the  bonding  charge  form  factors  depend  on 
the  amount  of  overlap  of  the  wavefunctions  centred  on  diff¬ 
erent  sites.  As  the  d  electrons,  unlike  tha  s  electrons, 
are  very  localized,  the  only  overlaps  that  would  contribute 
significantly  to  the  form  factors  are  those  arising  from 
the  hybridized  s-d  states. 


(3.22) 


(3.23) 
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Since  the  next  nearest  neighbour  is  not  much  further 
than  the  nearest  neighbour,  any  analysis  of  the  scattering 
factors  of  a  b.c.c.  transition  metal  must  include  at  least 
the  first  two  superlattice  form  factors.  Then,  writing 
equation  (3.21)  explicitly,  we  have 

S(G)  =  F  (G)  l  exp(-iG.Ri)  +  F^G)  £  exp(-iG.Si> 

1  x 

+  F2(G)J  exp(-iG.Ti)  (3.24) 

l 

where  is  the  position  vector  cf  each  basis  atom  and 
and  sre  the  position  vectors  of  the  scattering  centres  of 
the  first  two  super lattices  in  the  unit  cell.  Thus  the 
scattering  factors  at  the  first  five  reflections  are 


S( 

100) 

=  2Fo  - 

2F2 

S( 

200) 

=  2Fo  - 

8F1  + 

6F2 

S( 

211) 

i 

0 

u. 

CM 

II 

2F2 

S( 

220) 

=  2Fo  + 

8Fi  + 

6F2 

S( 

310) 

=  2Fo  - 

2F2 

It  is  clear  that  2FQ  will  be  the  scattering  factor  of  the 
free  atom  and  that  F^  and  Fj  will  .be  the  fourier  transforms 
of  the  products  of  hybrid  orbitals  centred  on  neighbouring 
sites.  Assuming  Slater  type  orbitals  for  the  s  and  d  func¬ 
tions,  we  see  that  the  bond  charges  can  be  represented  by 
exponen Lially  decaying  or  Gaussian  type  functions  centred 
at  the  bond  centre.  The  fourier  transforms  of  such  func¬ 
tions  will  also  be  decaying  functions,  with  a  maximum  at 
G  =  0,  and  will  be  positive,  at  least  for  small  G. 

The  scattering  factors  of  Fe  have  been  measured  by 
Eatterman  et  al.  (1960)  and  by  Paakari  and  Suortti  (1967). 
Their  results  together  with  the  Hartree-Fock  free  atom  cal¬ 
culations  for  the  first  five  reflections  are  shown  in  table 
3.1.  It  now  appears  that  Batterman  et  al.  did  not  take  pro¬ 
per  account  of  extinction  effects  in  their  analysis  of  the 
scattering  intensities.  The  deviation  of  their  results 
from  the  free  atom  calculations  are,  however,  confirmed 


TABLE  3.1 


Absolute  experimental  and  theoretical  atomic  scatt¬ 
ering  for  Fe. 

hkl  sin  8/ A 

110  0.247 

200  0.349 

211  0.427 

220  0.493 

310  0.551 

(a)  Paakari  and  Suortti  (1967) 


f  .  (a)  f.,  (free  atom) 

_ expt  •  tn 


18.19 

18.51 

15. 1910.08 

15.27 

13.01l0.08 

13.13 

11. 6010. 05 

11.61 

10. 4710. 05 

10.49 
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although  to  a  lesser  degree,  by  Paakari  and  Suortti.  The 
latter  also  found  the  deviation  to  be  negligible  beyond 
the  fifth  reflection. 

From  equation  (3.25),  we  see  that  the  deviation  bet¬ 
ween  the  experimental  and  free  atom  calculated  values  can 
be  accounted  for  solely  by  the  second  superlattice  bond 
charge.  We  have,  therefore,  plotted  these,  from  the  meas¬ 
urements  of  Paakari  and  Suortti,  against  sin  9/X  =  G/4*  in 
figure  3.1.  Although  the  errors  are  quite  large,  we  can 
easily  fit  a  smooth  curve  through  the  data  as  would  be  req¬ 
uired  if  the  bond  charge  model  were  a  good  description  of 
the  electronic  states  in  Fe.  Extrapolating  the  contribut¬ 
ions  of  Fq  and  F^  S[20(f]  and  Sp20j  from  table  3.1  and  fig¬ 
ure  3.1,  and  substituting  them  into  equation  (3.9),  we  find 
that  F^(k)  is  more  localised  and  smaller  than  FjOO.  This 
means  that  the  charge  distribution  around  the  first  neigh¬ 
bour  bond  centre  is  less  localised  and  less  important  than 
that  around  the  next  nearest  neighbour  bond  centre.  This 
is  consistent  with  the  fact  that  the  real  space  volume  is 
greater  in  the  jjLOCf]  direction  interstitial  site  than  in  the 
[ill]  one. 

S  3 . 4  Pi s  cussion 

Although  we  have  only  considered  Fe,  the  scattering 
factors  of  the  other  b.c.c.  transition  metals  also  show  a 
marked  discrepancy  between  the  experimental  and  free  atom 
values.  Cooper’s  (1964)  experiment  on  Cr  also  gave  large 
discrepancies  but  like  that  of  Batterman  et  al.  on  Fe  is 
probably  in  error  due  to  incorrect  extinction  corrections. 
Recent  experiments  of  V  (Mazzone  and  Diana  1978,  Korhonen, 
1971)  have  shown  the  discrepancy  at  the  first  reflection  is 
comparable  to  that  of  Fe,  so  we  can  have  confidence  in  res¬ 
ults  of  Paakari  and  Suortti. 

Weiss  (1978)  attempted  to  resolve  this  discrepancy  in 
the  b.c.c.  metals  by  contracting  the  3d  triply  degenerate 
radial  wavef unctions ,  which  point  to  the  nearest  neighbours 
and  expanding  the  3d  doubly  degenerate  radial  wavef unctions , 


which  point  to  the  second  nearest  neighbours. 

However,  he  had  only  partial  success  as  some  discrep¬ 
ancy  between  the  calculated  and  experimental  values  of  the 
[lltf]  scattering  factor  of  V  remains.  According  to  our 
formulation  of  the  scattering  factors,  his  model  only  adds 
angularity  to  the  one  centre  terms  and  does  not  consider 
the  two  centre  ones.  The  remaining  discrepancy  is,  there¬ 
fore,  due  to  these  two  centre  contributions. 

This  discrepancy  between  experimental  and  free  atom 
scattering  factors  is  also  found  in  the  f.c.c.  transition 
metals,  although  it  is  smaller  and  is  finite  only  for  small 
G.  By  a  similar  analysis  to  that  performed  in  §3.3,  we  can 
show  these  discrepancies  are  also  attributable  to  the  bond 
charges.  The  localisation  of  the  form  factors  of  the  bond¬ 
ing  charges  in  the  f.c.c.  metals  is  indicative  of  delocal¬ 
ised  bonding  electrons. 

The  present  approach  will  not  be  favourable  when  the 
WF's  are  long  range,  because  then  one  needs  to  consider  the 
overlaps  between  many  centres.  It  will  be  best  for  crystals 
in  which  even  the  valence  electrons  can  be  described  by  the 
tight  binding  approximation.  In  the  transition  metals,  if 
we  u3  3'u;r.c  only  overlaps  between  the  hybridised  orbitals,  we 
need  only  consider  overlaps  upto  the  second  nearest  neigh¬ 
bour.  Beyond  that,  the  size  of  the  overlaps  are  negligible. 

The  experimentally  determined  X-ray  scattering  factors 
then  tell  us  that  the  bonding  charges  are  included  in  the 
WF  description  and,  thus,  there  is  no  conflict  between  the 
"band"  picture  and  that  based  on  chemical  bonding.  Further, 
it  appears  that  bonding  is  less  prevalent  in  the  f.c.c. 
metals  than  in  the  b.c.c.  metals,  justifying  our  approxi¬ 
mation  in  Chapter  2.  From  our  analysis  of  the  scattering 
factors  we  find  that  second  neighbour  bonding  is  more  impor¬ 
tant  than  nearest  neighbour  bonding  in  the  b.c.c.  trans¬ 
ition  metals,  resulting  in  a  localisation  of  electrons  in 
the  [_10C[]  interstitial  site.  This  is  consistent  with  Weiss 
who  extended  the  radial  distribution  along  the  next  nearest 
neighbour  direction  and  contracted  it  along  the  nearest 


neighbour  direction. 


The  electronic  charge  distribution  in  the  b.c.c.  tran¬ 
sition  metals  can,  therefore,  approximately  be  written  in  the 
form: 

'  l  4ion  (£-£i>  *  l  4bond  1  <£'Si) 

1  1 

*  I  4bond  2  <3’26> 

where  R^,  and  T\  have  the  same  meaning  as  in  equation 

(3.24)  and  A.  ,  A,  ,  and  A.  ,  0  refer  to  the  ionic, 
ion’  bond  1’  bond  2 

nearest  and  next  nearest  neighbour  bond  charges  respectiv¬ 
ely.  We  have  used  this  model  to  calculate  the  phonon  freq¬ 
uencies  of  Fe  in  Chapter  5. 
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CHAPTER  4 


X-ray  Scattering  in  Be 


54.1  Introduction 


In  the  last  chapter,  we  saw  that  the  WF  description 
of  the  valence  electrons  was  very  useful  in  understanding 
the  measured  X-ray  scattering  factors  in  crystals.  Since 
the  scattering  factor,  in  this  formalism,  is  written  as  a 
sum  of  fourier  transforms  of  the  overlaps  of  the  WF's  cen¬ 
tred  on  different  lattice  sites,  the  expression  is  quickly 
convergent  only  if  the  WF's  are  short-range.  Hence,  this 
approach  is  best  for  those  crystals  in  which  the  valence 
electrons  are  best  described  by  tight-binding  wavefunctions 
However,  the  big  problem  is  then,  in  spite  of  the  value  of 
the  formalism,  the  practical  calculation  of  the  WF's. 


To  make  a  quantitative  test  of  our  expression  for  the 
scattering  factors  we  have  chosen  h.c.p.  Be.  This  choice 
is  based  on  the  following  considerations: 

(i)  the  valence  electrons  of  Be  as  calculated  by 
Donovan  (1952),  appear  to  be  tight-binding  in  nature. 

Brown  (1972)  in  her  analysis  of  the  scattering  factors  also 
arrived  at  the  same  conclusion. 

(ii)  the  tightbinding  nature  of  the  valence  electrons 
and  the  small  core  render  it  possible  to  obtain  useful 
approximations  of  the  WF's. 

(iii)  Inoue  and  Yamashita  (1973)  in  an  APW  calculation 
on  Be  estimated  that  both  the  lattice  parameters,  c  and  a, 
in  the  h.c.p.  lattice  would  have  to  be  increased  by  over 
160%  before  it  became  insulating.  This  implies  an  apprec¬ 
iable  overlap  between  the  2s  and  2p  bands  in  the  metallic 
phase  at  atmospheric  pressure,  which  will  result  in  substan 
tial  hybridization  between  these  bands. 

(iv)  Hybridization  leads  to  greater  overlap  and  thus, 
corrections  to  the  full  BZ  scattering  factor  are  not  neglig 
ible . 

(v)  the  X-ray  results  of  Brown  (1972)  indicate  subst¬ 
antial  angularity  in  the  charge  distribution. 
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Altmann  et  al.  (1957)  proposed  that  the  geometrical 
crystal  structures  of  the  d-band  metals  are  related  to  the 
formation  of  hybrids  consistent  with  these  structures. 

Yang  and  Coppens  (1978)  extended  this  idea  and  suggested 
that  hybrids  which  have  the  crystal  symmetry  are  also  pres¬ 
ent  in  h.c.p.  Be.  By  applying  group  theoretical  techniques 

o 

it  can  be  shown  that  in  the  h.c.p*  structure,  the  sp 
(s,  Px>  Py)  trigonal  hybrids  account  for  next  nearest  neigh¬ 
bour  planar  bonding  and  the  sp  (s,  p  )  linear  hybrids  for 
nearest  neighbour  interplanar  bonding.  The  charge  density 
distribution  maps  are  also  found  to  be  consistent  with  this 
bonding  scheme  (Yang  and  Coppens,  1978). 

In  h.c.p.  Be  the  c/a  ratio  is  much  less  than  that  of 
the  ideal  h.c.p.  lattice.  With  the  above  bonding  scheme, 
we  see  that  this  indicates  stronger  interplanar  bonding 
than  that  within  the  planar  layers,  which  is  consistent 
with  the  results  of  Brown  (1972). 


In  our  attempt  to  calculate  the  scattering  factors  of 
Be,  we  have  made  the  following  assumptions  and  observations: 

a)  The  valence  band  is  a  mixed  s-p  band. 

2 

b)  sp  bonding  is  so  much  Jess  prevalent  than  sp 
bonding  that  it  can  be  neglected.  This  is  very  helpful  in 
the  construction  of  the  WF’s. 

c)  The  linear  sp  hybrids  can  be  shown  to  form  the 
orthonormal  set 


1  (s 

/? 


+  pz) 


and 


1  (s 

n 


pz 


) 


which  we  assume  to  be  equally  populated.  s  and  p  refer  to 
the  2s  orbital  and  z  component  of  the  2p  orbital  respectiv¬ 
ely,  where  z  is  taken  to  be  perpendicular  to  the  hexagonal 
plane . 

d)  The  px  and  Py  bands  are  unoccupied.  As  the  number 
of  Bloch  waves  is  equal  the  number  of  WF  bands,  the  sp  band 
has,  therefore,  four  WF's  of  which  only  two  are  occupied. 

In  the  next  section  we  define  WF’s  in  a  complex  latt¬ 
ice  and  then  modify  the  expression  for  S(G)  for  a  h.c.p. 
lattice  appropriately.  By  making  use  of  the  above  assump- 
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tions,  WF's  are  then  constructed  for  Be  in  §4.3. 

In  §4.4  we  report  the  calculated  scattering  factors 
of  Be  and  compare  them  with  the  experimental  results.  The 
validity  of  certain  approximations  that  we  have  made  and 
the  actual  results  are  discussed  in  §4.5. 


§4.2  Wannier  Functions  and  Scattering  Factors  in  Complex 
Lattices 


Complex  lattices  are  those  lattices  in  which  the  unit 
cell  contains  more  than  one  atom.  Then,  by  analogy,  we  def¬ 
ine 

,  (r)  =y7-j~y  l  exp  (ifc.R  )  a  (r-R  )  (4.1) 

ok  -  / ( Ns )  L  *  —  — p  a  —  — p 

—  p 

and 

aa(£_V  =  /dt)  l  exp  *«k  (4*2) 

These  functions  now  satisfy  all  the  usual  orthonormality 
conditions  which  have  to  be  obeyed  by  n^och  functions  and 
WF's. 


In  Be  there  are  two  atoms  per  unit  cell  and  so  we  can 

consider  the  cell  to  be  a  mesh  of  two  sublattices.  If  we 

denote  the  position  vectors  of  the  atoms  of  one  sublattice 

by  Rn,  t^ose  of  the  other  can  be  denoted  by  Rn  +  d.  Because 

both  atoms  are  related  through  a  centre  of  inversion,  it  is 

convenient  to  choose  local  axes  x,  y  and  z  on  the  second 

atom  inverted  with  respect  to  those  on  the  first.  Then,  if 

the  WF's  centred  on  the  atoms  of  one  sublattice  are  denoted 

by  a(r  -  R  ) ,  the  WF's  centred  on  the  atoms  of  the  other 

sublattice  will  be  a(R  +  d  -  r) 

— n  —  — 

Hence,  omitting  the  band  index,  from  equation  (4.1) 
the  Bloch  function  is  given  by 

=  *TnI7  I  {a(r-Rn)  +  a(Rn+d-r) exp( ik. d) }exp( ik.R^) 


(4.3) 


Then,  by  a  similar  analysis  to  that  carried  out  in  53.2,  we 
find  that  the  structure  factor  can  be  written  in  the  form 


S(k)  =  7  y  exp( ik . (R-R  ) )  [a  (r-R  )a  (r-R_)exp(  ik  . r)dr 
—  r  —  — m  — n  a  —  — n  e  —  — m  *  —  —  — 

a  B  kmn  1  p 

+  y  y  exp(ik.(R  -R  ))  [a  (R  +d-r)a„(R  +d-r)exp(ik.r)dr 

" 0  t  —  — m  — n  I  a  — n - g  — m - r - — 

a  B  kmn  1 

+  y  y  exp(ik.(R  -R  ))  a  (r-R  )a.(R  +d-r)exp( ik . r)dr 

L0  b _ e  —  — m— n  I  a  —  —m  B  — n - e - — 

a B  kmn  J 

y  y  exp( ik . (R-R  ) )  [a  (r-R  )a.(R  +d-r)exp(ik .r)dr 
L  b  r  —  — m  — n  I  a  —  —  n  B  —m  —  —  c  —  —  — 

aB  kmn  ;  p 


(4.4) 


where  a  and  B  refer  to  the  (s  +  p  )  and  (s  -  p  )  bands  res- 

z  z 

pectively . 


We  now  make  the  following  substitutions  in  equation 

(4.4) 

(i)  In  the  first  sum  r  r  +  R  +  R„ 

2 

(ii)  In  the  second  sum  r->-r  +  R  +  R  +2d 

—  —  — m  — n 

2 

(iii)  In  the  third  and  fourth  sums  r-*r  +  R„+R+d 

2 

Then,  it  follows  that  the  scattering  factors  are  zero  except 
at  the  reciprocal  lattice  vectors,  where  they  are  given  by 
the  expression 


S  (G ) 


y  3;5l(kFRn)  y  exp(-iG.(Rn  +  f))  I_a(R_,G) 

aBn  “O -  f  “  T  ~  -n  - 

r  n 


+  y  3jl(kFRnf)  y  exp(-iG.(^n  +  f))  i  (RJ^,G) 
aBn  kFRn'  f  aB  n 

(4.5) 


where 
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and  f  are  the  position  vectors  of  the  atoms  in  the  sub¬ 
lattice.  For  Be ,  f  =0  or  d. 

The  overlap  integrals  I(Rnj£)  and  I(Rn',G)  are  def¬ 
ined  by  the  equations 

f  R  R 

Ia3^— n’— ^  =  j  aa^— ^  a^(r  +  ?)  exP  ^iG.r)dr  (4.6) 


which  is  the  fourier  transform  of  the  overlap  between  fun¬ 
ctions  centred  on  atoms  of  the  same  sublattice,  and 

r  R  *  -R  ' 

IagCRn’.G)  =  I  aQ ( r-  ~~)  a^(—~ - r)  exp(iG.r)dr  (4.7) 

which  is  the  fourier  transform  of  the  overlap  between  fun¬ 
ctions  centred  on  atoms  of  different  sublattices. 


S 4 . 3  Construction  of  the  Wannier  Functions 

As  mentioned  in  the  Introduction  to  this  Chapter,  we 
can  construct  approximate  WF's  for  Be  because  of  the  tight 
binding  nature  of  the  valence  electrons.  Comparing  the 
tight  binding  Bloch  function  to  that  given  by  equation  (4.3) 
we  see  that  a  first  order  approximation  of  the  WF's  are  in 
fact  the  atomic  orbitals.  However,  due  to  the  orthogonality 
property  of  the  WF’s,  the  atomic  orbitals  have  to  be  made 
orthogonal  to  each  other. 

In  Be  we  assume  only  two  occupied  bands.  If  we  den¬ 
ote  the  (s  +  pz)  atomic  orbital  centred  on  Rn  by  ^^(r-Rn) 
and  the  corresponding  (s  -  pz)  atomic  orbital  by  ^2^— n^ 
then 


.  (m,n) 

— 


ip.(r-R  )  *|».(r-R  )  dr 
J  1  —  — m  j  —  — n  — 


(4.8) 


is  the  (i  +  m)  by  (j  +  n)  overlap  matrix  between  these  fun¬ 
ctions  centred  on  different  lattice  sites. 
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If  the  ip's  are  WF's,  we  require 


A . . (m,n) 


6  .  .  6 
l  j  mn 


(4.9) 


To  satisfy  equation  (4.9),  we  need  to  orthogonalise 
the  tp's.  We  transform  the  set  {4/(r-Rm)}^  to  an  orthog¬ 
onal  set  (Wr-R^J^  by  the  symmetric  orthogonalization 
procedure  introduced  by  Ltiwdin  (1956). 


It  can  be  shown  that  the  complete  orthogonal  sets  can 
be  written  in  terms  of  the  non-orthogonal  atomic  functions, 
nemely 

W1<£-V  *  iii1/2  Vn-V 

and 

W0 (r-R  ) 

2  —  — m 


where  A12*  —21  ancJ  —22  are  rnatrices  making  up  the 

block  overlap  matrix  A^(m,n)  =  A  .  From  equations  (4.8)  and 
and  (4.9)  we  see  that  for  the  WF's  to  be  exact  the  A  mat¬ 
rices  are  infinite.  This  is  equivalent  to  saying  that  the 
WF  on  one  lattice  site  is  orthogonal  to  WF's  on  all  other 
lattice  sites.  This  is  of  course  not  practical  in  any  cal¬ 
culation.  The  accuracy  of  the  constructed  WF's  is,  there¬ 
fore,  dependent  on: 

(i)  the  order  to  which  the  a.,  matrices  are  taken; 

(ii)  the  order  to  which  the  inverses  and  products 
in  equations  (4.10)  and  (4.11)  are  calculated;  and 

(iii)  the  size  of  the  elements  of  the  a  matrices. 

Hence,  equations  (4.10),  (4.11)  and  (4.8)  provide 
the  basis  of  an  approximate  calculation  of  the  WF's. 


—22  (»2(^m)  ~  ill  il2*l(£-gm)) 

(1  "  ill  i,22  ^21  -12) 


(4.11) 


14.4  Calculation  of  the  X-ray  Scattering  Factors  in  Be 

Be  has  a  h.c.p.  crystal  structure  with  lattice  para¬ 
meters  a  =  2.28588  and  c  =  3.58438.  In  the  calculation  of 


of  the  WF's  we  require  the  atomic  2s  and  2p  orbitals,  which, 
for  our  calculations,  we  have  taken  to  be  the  orthogonal- 
ized  Slater-type  analytic  functions  given  by  Duncanson  and 
Coulson  (1944).  The  screening  parameter  of  the  2p  orbital 
was  extrapolated  from  that  given  for  the  heavier  elements. 

The  normalised  hybrid  wavefunctions  then  are 

lp^(r)  =  ^  (Ar  exp  (-yr)  -  B  exp  C-vr)  +  Cz  exp  (-Xr)) 

(4.12) 

A 

=/2  ^Ar  exP  ^”4r)  -  B  exp  (-vr)  -  Cz  exp  (-Xr)) 

(4.13) 

where  the  constants  are  given  in  table  4.1. 

Since  the  only  hybrids  we  allow  in  this  system  are 

the  (s  ±  p_)  hybrids,  we  need  only  consider  the  overlaps 
z 

between  the  wavefunctions  centred  on  one  atom  and  its  neigh¬ 
bours  up  to  and  including  the  third  shell.  This  corresponds 
to  taking  its  first  fourteen  neighbours. 


For  the  (s  ±  p  )  hybrid  orbitals  the  overlap  matrices 
z 

defined  in  §4.3  can  be  written  in  terms  of  the  following 
six  overlaps: 

°1  =  j  ^1^—^  =  {  ^2^—^  ^2^— 5^d— 

°2  =  }  ♦i^  ^(r-Rg )dr  =  J  if>2(r)  8^d— 

a 3  =  |  ^(r)  ^  (r-I^dr  =  J  ^(r)  'f'2^-_^.S>d- 

04  =  }  ^  (r)  ip 2  (r-Rj  Mr  =  J  4>2  (r)  >f)^(r-Rs)dr 

a5  =  |  i|> ^(r)  i)/2 (r-Rg )dr  =  J  i|>2(r)  ^1(r-Rg)dr 

0g  =  |  ^(r)  'J>2<r-^4)dr  =  J  *2(r)  ^(r-R^dr 


where  Rn  are  the  position  vectors  of  the  atoms  labelled  by 
n  relative  to  the  ’zero’  atom  labelled  1.  From  the  spatial 
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TABLE  4.1 


Parameters  of  2s  and  2p  wavefunctions  of  Be. 


Parameter 

Value 

a  (r5/2> 

0.837 

B  <r3/2> 

1.788 

c  <r5/2> 

1.128 

o—  1 

V  (A  ) 

1.851 

o-l 

v  (A  ) 

6.200 

o-l 

X  (A  ) 

1.701 
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extent  of  the  hybrid  orbitals,  it  is  clear  that  all  other 
overlaps  will  be  much  smaller  and  hence,  can  be  neglected. 

Using  the  functions  given  in  equations  (4.12)  and 
(4.13),  we  have  calculated  these  overlaps  and  they  are 
presented  in  table  4.2.  From  their  definition,  it  follows 
that 


4n  ~  2 

and 

£21  =  £.12 

where  the  dagger  t  denotes  the  hermitian  conjugate.  For 
the  system  of  fifteen  atoms  considered,  we  can  express  the 
^  matrices  as 

A.  1  =  1,  +  S  (4.15) 

and 


where  S  and  T  depend  on  the  o^(i  =  1  to  6). 

Due  to  the  complexity  of  the  S_  and  T  matrices,  we  have 
worked  only  to  order  ST  and  TT+  in  the  expansion  for  the 
approximate  WF's  from  equations  (4.10)  and  (4.11).  The 
validity  of  this  approximation  is  discussed  in  the  next 
section. 


It  may  be  easily  shown  that 
W1(r)  =  ^(r)  -  —■  (C1(£)  +  ^(r))  -  ~  C1(r) 

~T-  <^1<I-^14)  +  ♦i<S-Ris>> 


(4.17) 


and 


W2(r)  = 


r  >60  - 


1  +  3 o*  +  6og  +  Ogj  ^(r) 


+  l  (2 o4o5  +  o4o6  -  a1)(C2(r)  +  n2(r)) 


+  }  (oj  +  2ag  -  a2)  ?2(r)-  ~  + 

+  Oc^c^  +  6o2o5  +  a3°g)  ^(r)  +  ■f(2oi°5  +03°4^1^ 
+  (|(o1ag  +  2 a1as  +  2o2<j4)  -  o4)  ^(r) 


+  (|(o1o4  +  2o2o5)  -  a5>  ^(r) 
+(|  o1o4  -  Og)  *1(r-R15) 


(4.18) 


where 


€•  =  J  i>  ■  (r-R) 


n .  =  j  ( r-R  ) 

l  i - n 

n  =  5 


C  =  I  .  (r-R  ) 

%  ri - n 

n=  8 

Since  we  have  only  considered  overlaps  up  to  and  inc¬ 
luding  the  third  shell  of  neighbours,  we  need  consider  only 
the  first  three  superlattice  contributions  to  the  scatter¬ 
ing  factor  in  equation  (4.5).  The  scattering  factor  at  the 
reciprocal  lattice  vectors,  defined  by  Miller  indices  ]hRl]  , 
is  then 


S[Ncl]  =  ls0B*f]  I 

*  a  6 


+  33l(kFd)  S1[hkl]  I  Ia6(d’hkl) 


+  3jl(kFa)  S  [hkl]  l  I  (a,hki: 


+  3jl(kFC>  S  [hkl]  l  I  B(c",hki: 

k  c  a6 

F 


(4.19) 


where  S  [hkl]  is  the  geometrical  structure  factor  of  the 
t  ^  -  _ 

n  superlattice,  IQg(x,hkl)  are  the  average  overlap  inte¬ 
grals  defined  by  equations  (4.6)  and  (4.7)  over  and  R  ' 
with  moduli  x,  and  d,  a  and  c  are  the  first,  second  and 
third  neighbour  shell  distances  respectively. 

The  X-ray  scattering  factors  of  Be  as  measured  by 
Brown  show  considerable  disagreement  with  the  free  atom  cal¬ 
culations  (Weiss  1968).  The  discrepancies  between  the  two 
sets  are  largest  at  the  first  five  observed  reflections. 

We  have,  therefore,  calculated  S(G)  from  equation  (4.19)  at 
these  reciprocal  lattice  vectors. 

The  WF's  were  constructed  by  the  method  described  in 
54.3.  These  were  then  used  in  equations  (4.6)  and  (4.7)  to 
evaluate  the  average  overlap  integrals  of  equation.  (4.19). 
Implicit  in  equation  (4.19)  is  the  fact  that  we  have  used  a 
spherical  Fermi  surface  approximation.  The  validity  of  this 
approximation  is  also  discussed  in  the  next  section. 

The  calculated  scattering  factors  are  tabulated  in 
table  4.3,  together  with  the  experimentally  determined 
values  of  Brown.  Also  shown  are  those  calculated  by  Brown 
using  a  tight  binding  approach,  and  the  free  atom  core  and 
valence  electron  contributions. 


54. S  Discussion 

Before  discussing  the  results  of  our  calculations  which 
are  presented  in  table  4.3,  we  make  a  few  comments  on  the 
two  major  approximations,  other  than  the  assumptions  of  the 
model,  that  we  have  used  in  the  calculations. 
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The  WF's  used  to  calculate  the  structure  factors  were  ob¬ 
tained  by: 

(i)  orthogonalizing  the  hybrid  orbitals  on  one  atom 
with  only  the  first  three  shells  of  neighbours; 

(ii)  truncating  the  resulting  expressions  to  order 
ST  and  T+T  of  the  overlap  matrices. 

Both  these  approximations  lead  to  an  error  in  the  long 
range  part  of  the  WF.  The  first  term  in  equation  (4.19)  is 
just  the  fourier  transform  of  the  square  of  the  WF .  Any 
errors  in  the  long  range  part  of  the  WF  will,  therefore, 
only  show  itself  at  small  G,  the  largest  error  being  at 
G  =  0.  We  estimate  that  even  at  the  first  reflection,  G  is 
sufficiently  large  to  minimize  the  errors  which  are  of 
order  1/Gr  where  r  is  the  interatomic  spacing.  The  other 
terms  in  equation  (4.19)  are  essentially  the  fourier  trans¬ 
forms  of  symmetric  functions  centred  away  from  the  origin. 
Hence,  once  again  the  errors  in  the  WF’s  will  not  lead  to 
any  significant  errors  in  the  scattering  factors  away  from 
G  =  0.  So,  although  we  can  make  no  claims  on  the  accuracy 
of  the  calculated  WF's  the  approximations  used  in  obtaining 
them  appear  to  useful  in  the  calculation  of  the  X-ray  scat¬ 
tering  factors. 

The  other  majcr  approximation  we  have  made  is  that  of 
a  spherical  Fermi  surface.  In  Be,  this  is  not  a  good  approx 
imation.  As  we  saw  in  Chapter  3,  the  first  term  in  equation 
(4.19)  is  that  which  would  be  obtained  for  an  insulator,  i.e 
one  with  a  filled  BZ  and  where  the  angularity  in  the  charge 
distribution  is  incorporated  into  the  WF.  Be  is  a  very 
good  example  of  a  compensated  metal  and  so  the  other  terms 
in  equation  (4.19'  are  just  those  contributions  arising 
from  the  holes  in  the  first  BZ ,  and  from  the  electrons  out¬ 
side  this  zone.  The  approximation  we  have  made,  therefore 
entails  describing  these  electron  and  hole  states  by  hybrid 
wavefunctions  with  |kj  =  k^,.  In  doing  this,  we  are  of 
course  giving  it  more  weight  than  is  due  because,  in  addit¬ 
ion  to  the  hybrid  states,  there  exists  at  this  wave-vector 
free  electron  like  states  which  will  contribute  differently 
to  the  scattering.  We  have  followed  the  above  approach  and 
not  used  plane  waves  for  the  following  reasons: 

(i)  plane  waves  (including  standing  waves)  have 
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fixed  values  of  k  and  are  hence  spatially  spread  over  the 
whole  volume  in  real  space.  This  does  not  take  account  of 
the  fact  that  these  states  must  be  excluded  from  the  core 
regions . 

(ii)  the  overlap  electrons  with  energy  near  will 
not  be  properly  represented. 

(iii)  to  find  its  correct  contribution  to  the  scattering 
factors,  we  would  first  have  to  calculate  the  contributions 
from  the  zones  into  which  the  electrons  have  spilled  over. 

In  Be,  this  would  include  the  second  and  third  BZ. 

2 

From  table  4.3.  we  see  that  the  free  atom  Is  core 
constantly  overestimates  the  scattering  factor  at  all  the 
reflections.  With  the  exception  of  the  first  reflection, 
we  find  that  the  introduction  of  angularity,  through  the 
hybridised  WF's  into  the  core  charge  density  and  the  incl¬ 
usion  of  the  overlap  terms  lead  to  smaller  scattering  factors 
than  that  predicted  from  free  atom  calculations.  Thus,  the 
general  trend  is  to  bring  the  theoretically  calculated  scat¬ 
tering  factors  more  into  agreement  with  experiment.  Our  res 
ults  are  also  in  good  agreement  with  the  tight-binding  cal¬ 
culations  of  Brown.  This  vindicates  our  model  with  the  ass¬ 
umption  of  only  sp  linear  hybridization  and  shows  that  most 
of  the  localised  charge  is  interplanar. 

Though  four  of  the  five  calculated  scattering  factors 
show  excellent  agreement  with  experiment,  the  one  serious 
disagreement  is  at  the  first  reflection.  This  is  not  a  res¬ 
ult  of  the  above  mentioned  approximations. 

If  we  consider  a  charged  'blob'  model  similar  to  that 
proposed  in  Chapter  3,  we  arrive  at  the  same  conclusions. 

In  such  a  model,  placing  charges  at  the  mid-points  between 
an  atom  and  its  first  three  shells  of  neighbours,  we  find 
that  the  scattering  factors  at  the  first  five  reflections 
are  given  by 
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S(100)  =  f0  +  2f1  -  f2  +  fg 

S( 002 )  =  2 ( f Q  -  3f±  +  3f 2  +  f  ) 

S(101)  =  3(fn  -  f,  -  f,) 

0  1  6  (4.20) 

S ( 10 2 )  =  f0  -  2f1  -  f2  +  fg 

S(110)  =  2 ( f g  -  f1  -  f2  +  fg) 

where  f  are  the  form  factors  of  the  ’blobs'  of  the 
n 

superlattice. 

Since  these  form  factors  are  fourier  transforms  of 

exponential  or  Gaussian  functions,  they  will  be  positive, 

at  least  for  small  G.  Also,  as  the  linear  sp  bonding  is 

o 

more  prevalent  than  the  sp  planar  bonding,  the  first  super¬ 
lattice  form  factor  will  be  greater  than  the  second  super¬ 
lattice  form  factor.  Hence,  from  equation  (4.20)  all  the 
scattering  factors,  except  at  the  first  reflection,  will 
be  less  than  the  free  atom  calculated  ones  which  is  in  agree¬ 
ment  with  our  calculations.  The  scattering  factor  at  Q.0cG 
will,  however,  be  larger  than  the  free  atom  result,  again 
in  agreement  with  our  calculations.  Hence,  the  discrepancy 
at  [lOCf]  is  not  a  result  of  the  approximations  made  in  the 
calculation,  but  something  more  subtle.  We  have  not  been 
able  to  f-'nd  a  satisfactory  explanation  but  suggest  it  might 
be  due  to  the  fact  that  we  did  not  take  account  of  the  plane 
wave  like  states  in  the  form  of  othogonalized  plane  waves. 
This  would  give  rise  to  a  form  that  is  finite  only  at  small 
G  (e.g.  at  the  first  reflection)  and  almost  zero  at  large  G. 

The  good  agreement  of  our  calculations  with  the  exper¬ 
imental  scattering  factors  at  G,  other  than  at  the  first 
reflection  justifies  the  model  that  we  have  used  and  the 
approximations  that  were  made.  Further,  it  shows  that  the 
WF  description  of  the  hybrids  is  a  useful  one,  and  also 
that  the  bond  charge  model  is  quite  realistic. 


CHAPTER  5 


The  interatomic  potential  in  the  B.C.C.  transition  metals 
§5.1  Tntroduction 

By  investigating  the  experimentally  measured  scatt¬ 
ering  factors  of  the  transition  metals,  we  saw  that  a  use¬ 
ful  representation  of  the  angularity  in  the  charge  distrib¬ 
ution  could  be  obtained  by  the  introduction  of  bond  charges. 

A  description  of  the  outer  electrons  in  terms  of  WF's  leads 
to  the  concept  of  bonding  charges  situated  at  interstial 
sites.  By  linear  response  theory,  these  give  rise  to  angu¬ 
larity  in  the  force  fields.  Thus,  as  angular  forces  are 
necessarily  required  in  a  description  of  the  transition 
metals,  any  model  must  in  some  way  incorporate  Trhese  bonding 
charges . 

From  our  discussion  in  Chapter  3  it  is  clear  that  the 
angularity  in  the  charge  distribution  and  hence,  in  the 
force  fields  is  greater  in  the  b.c.c.  than  f.c.c.  transition 
metals.  We  have,  therefore,  made  a  model  for  the  b.c.c. 
transition  metals  from  which  the  phonon  frequencies  and 
elastic  properties  may  be  calculated. 

Iii  our  model,  in  addition  to  the  ions  situated  at  the 
lattice  ..ites  we  place  bond  electrons  at  interstitial  sites 
consilient  with  the  symmetry  of  the  crystal.  The  bonds  are 
constructed  in  the  manner  proposed  by  Altmann  et  al. (1957).  In 
this  model  the  charge  density  will  then  be 


p(r)=  y  A.(r-R  )  +  T  A,(r-r  ) 
—  L  l  —  — n  L  b  —  — v 

n  v 


(5.1) 


where  R  denotes  the  ionic  positions  and  r  the  bond  elec- 
— n  —v 

tron  positions.  Since  the  bond  electrons  are  situated  half¬ 
way  along  a  bond  connecting  two  atoms  ,  we  can  write 


R  +  R 

r  =  — n  — m 
—v  - 


(5.2) 


As  the  first  two  superlattice  form  factors  in  Fe  are  suff¬ 
icient  to  explain  the  X-ray  scattering  factor  results  equat 
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ion  (6.1)  can  be  reduced  to  equation  (3.26).  Hence,  the 
picture  of  the  b.c.c.  transition  metal  is  that  of  ions  sit¬ 
uated  at  the  lattice  sites  and  bond  electrons  situated  at 
the  first  two  superlattice  sites  all  immersed  in  an  uniform 
electron  gas.  The  ions  then  interact  with  each  other  dir¬ 
ectly  and  via  the  electron  gas. 

This  interaction  takes  the  same  form  as  for  the  sim¬ 
ple  metals.  However,  in  addition  to  this  interaction, 
there  are  also  the  interactions  of  the  bond  electrons  bet¬ 
ween  themselves  and  with  the  ions.  These  interactions  will 
also  be  screened  by  the  electron  gas.  However,  their  exact 
nature  has  not  been  investigated. 

If  we  denote  the  ion-ion,  bond  electron-ion  and  bond 
elec Iron-bond  electron  interactions  by  0^,  0?  and  0^  res¬ 
pectively,  the  total  energy  of  the  system  can  be  written  in 
the  form 


where  Lho  set  {ry}  are  completely  determined  by  the  set 
{ Rn }  through  equation  (5.2)  with  the  condition  that 


|R  -R  I  =  R 

where  Rc  is  the  first  or  second  neighbour  distance. 

In  writing  this  we  have  assumed  that  the  bond  charges 
are  spherically  symmetric  about  the  bond  centre  ensuring 
that  the  bond  electron-ion  and  bond  electron-bond  electron 
interactions  are  central. 


In  §5.2  we  derive  expressions  for  the  second  order 
elastic  constants  and  the  phonon  frequencies  in  cubic  metals 
with  metallic  binding. 

These  expressions  are  used  in  15.3  to  calculate  the 
second  order  elastic  constants  and  phonon  frequencies  of 


Fe.  We  have  chosen  Fe  as  an  example  because: 


(i)  the  X-ray  scattering  factors  show  that  its  charge 
distribution  is  quite  angular  and  can  be  written  in  the  form 
given  in  equation  (3.26). 

(ii)  a  liquid  pair  potential  has  been  calculated  by 
Waseda  and  Ohtani  (1975b).  This  provides  a  starting  point 
in  the  determination  of  0^(r)  in  equation  (5.3).  The 
liquid  potential  was  then  refined  and  the  effects  of  the 
bond  charges  on  the  phonon  examined.  Because  of  the  simil¬ 
arity  of  the  Cauchy  discrepancy  and  the  phonon  dispersion 
curves  with  that  of  K,  it  would  appear  that  the  forces  in 
Fe  are  almost  central.  However,  it  is  clear  from  the  argu¬ 
ments  presented  in  Chapter  3  that  this  is  not  the  case. 
Hence,  this  similarity  must  be  due  to  some  cancellation  bet¬ 
ween  the  effects  of  the  bond  charges  and  that  of  Ev>  In 
this  sense,  Fe  is  not  the  best  choice,  but  the  lack  of  a 
liquid  pair  potential  in  any  of  the  other  b.c.c.  transition 
metals  has  led  us  to  make  this  particular  choice. 

In  §5.4  we  discuss  the  usefulness  of  the  bond  charge 
model  and  ways  in  which  a  correct  description  of  the  force 
field-  an  be  obtained.  As  the  central  pair  potential  gives 
an  ircp'.i  lant  contribution  to  the  total  energy,  we  have 
investigated  its  effects  on  the  liquid  structure  factor. 


15.2  Expressions  for  the  Elastic  Constants  and  Phonon 
Frequencies 

In  the  following  derivations,  we  assume  that  the  bond 
charges  do  not  change  appreciably  under  an  infinitesimal 
deformation;  hence,  the  potential  functions  0j  an^  ^3 
remain  unaltered. 

a)  Elastic  Constants  for  a  Cubic  Crystal 

If  there  are  2J  bonds  pointing  out  from  one  atom,  then 
there  are  J  bond  electrons  associated  with  each  atom  with 
position  vectors  x°(a  =  relative  to  the  atom.  The 

bond  energy  per  unit  undeformed  volume  can  be  written  in  the 
form 
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u  =  2F"  l  K<rS)  +  2  I  f2(^S  + 

os'  a 

+  I  d3<  !rS  +  ia  -  l\  >]  +  ev(v) 

ft  T  ft 


(5.4) 


s  •  •  th 

where  R  are  the  position  vectors  of  the  s  atom  relative 

to  some  'zero'  atom  and  fi  is  the  undeformed  atomic  volume. 

o 

ev(V)  is  defined  by  equations  (2.14)  and  (2.38). 

Under  a  homogeneous  deformation,  we  let 

_  s  s  „  S  s 

R  -*•  r  =  R  +  u 


Then  adopting  a  similar  procedure  and  notation  to  that  used 
in  Chapter  2,  it  can  be  shown  that  the  first  and  second 
order  elastic  constants  in  a  cubic  crystal  are 

Cf  =  I  f*i(xS)2  +  l  ,'xS‘xS  +  l  l  'J,3xS ^xS]  +  =  0  (5.! 

s'-  a  aS?3  ' 

C11  =  lUi<xS>4  +  Jl*2(xS)2(xSa)2  +  *2{(xS)2-xSaxS)) 

"s'-  a 

+  4  £  (  i|j  3*(xS )  2(xS® )  2  +  ^3((xs)2-xSPxS))]  + 

(5.6) 


C12  =  i'K,(xS)2(yS)2  +  i  l  *2*VxSV" 

R  '  n 


,  r  s  s  SI  si  "  ' 

t  5  V  *3x  y  X  y  *  ♦  ev 

Of  P  ' 
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c44  =  I  f^1(xs)2(yS)2 

♦  £  l  2L''(Cxs>2<ysa>2  *  xVxSV°)  ♦  ^Hx^-tx^x*)) 

a 

*3!  f  2  |^3'(  (xs)  2(ysh2+  xVx6V6)  *  *2<(xs>2-2xeSx6>>) 


(5.8) 


The  Cauchy  discrepancy  given  by 
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+  ^'(2xsexs-(xs)2))  +  e"  +  2e 


(5.9) 


now  contains  both  structure  dependent  and  structure  indep¬ 
endent  quantities.  Equation  (5.9)  shows  explicitly  the 
angular  contribution  to  the  Cauchy  discrepancy. 


b)  Phonon  Frequencies 

The  phonon  frequencies  in  this  model  are  given  by 
equations  (2.26)  to  (2.28),  with  the  potential  energy  given 
by  equation  (5.3). 

After  performing  the  differentiations,  we  find  that  the 
contribution  from  the  structure  dependent  quantities  to  the 
dynamical  matrix,  (k)  is  given  by 


MD^(y 


(r(l'))  [1 


exp(ik.r(l'  )  f] 


+  *  \  ,  02ag 
1  s 

[r(s  1  )  =±2Ta} 

4-Jst  02aB 

[r(l,)=±2T°] 


iph- 

(PP- 


r(l')  1)0  -  exp( ik.r(  1 ')f] 
r(s)  |]  Q  -  exp( ik.r(l '  ))3 


*!£’,  «2..(|£(S'>^<1'-)-|)ll-  exp(iK.r(l')0 

1  s 

[r(s  '  )-r(  1  ’ )  =  ±21**]  (5.10) 


where 


Equation  (S.10)  together  with  the  contribution  from  Ev(V), 
given  in  §2.2  thus,  enables  us  to  calculate  the  phonon  freq¬ 
uencies  in  the  three  symmetry  directions  in  the  bond  charge 
model. 


§5.3  Pair  Potentials,  Phonon  Frequencies  and  Force  Fields 
in  Fe 

The  arguments  presented  in  Chapter  2  as  to  why  the 
central  pair  potential  in  the  solid  state  can  in  the  first 
instance  be  approximated  by  the  liquid  pair  potential, 
remain  valid  in  the  b.c.c.  transition  metals.  Waseda  and 
Ohtani  ( 1979a)  measured  the  liquid  X-ray  structure  factor  of 
Fe  and  extracted  a  pair  potential  by  the  Born-Green  (1946)  method. 
The  main  difference  between  the  solid  and  liquid  pair  pot¬ 
entials  will  be  the  way  in  which  it  falls  off  at  large  r. 

We  have  fitted  the  pair  potential  of  Waseda  and  Ohtani  to 
the  form  given  by  equation  (2.37).  The  procedure  was  the 
same  as  that  used  for  Ni.  The  hard  core  again  character¬ 
ised  by  a  repulsivity  of  sixteen  was  influenced  by  the  melt¬ 
ing  and  liquid  structure  factor  data. 

Our  interpretation  of  the  measurements  of  the  melting 
temperatures  under  pressure  in  §2.3,  lead  us  to  assume  the 
repulsivity  of  Fe  to  be  of  the  order  of  fifteen.  The  simil¬ 
arity  of  the  long  range  liquid  S(k)  of  Fe  and  Ni  reinforces 
our  choice  of  n.  The  major  difference  between  the  solid  and 
liquid  potentials  is  in  the  mean  free  path,X  .  In  the  liquid 
state,  due  to  the  disorder,  X  is  about  five  times  greater 
than  in  the  solid.  Since  the  nearest  neighbour  spacing  in 
Fe  is  almost  the  same  as  in  Ni ,  we  expect  the  inverse  mean 
free  paths  in  the  solid  and  liquid  to  be  almost  the  same 
namely  0.05&  *  and  0.20$  ^ . 

Using  the  liquid  potential,  but  with  as  the 

inverse  mean  free  path,  we  have  calculated  the  phonon  freq¬ 
uencies  in  b.c.c.  Fe.  The  parameters  in  the  potential  were 
adjusted  until  the  best  fit  to  the  experimentally  determined 
phonon  frequencies,  which  is  shown  in  figure  5.1,  was  obtained. 
The  resulting  central  ion-ion  pair  potential  was  given  by 
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01(r)  =  +  ~  cos  ( 2 . 36r  +  1.41)  exp  C-0.05r)eV 

r  r  (5.11) 

and  is  shown  in  figure  5.2.  In  calculating  the  phonon  freq- 

o 

uencies  ,  we  have  truncated  the  potential  at  7A. 

The  lack  of  knowledge  about  the  bond  electron-bond 
electron  and  bond  electron-ion  interactions  requires  us  to 
take  an  empirical  approach  in  incorporating  the  bond  charges 
into  the  formulation  of  the  total  energy. 

Since,  in  the  calculation  of  the  elastic  constants  and 
of  the  phonon  frequencies,  we  require  only  the  first  and 
second  derivatives  of  the  interacting  potentials  at  the 
lattice  and  bond  centre  sites,  we  have  treated  them  as  para¬ 
meters  . 


If  the  bond  charges  are  small  compared  to  the  ionic 
charges,  it  follows  that  0j  will  be  larger  than  0^  in  mag¬ 
nitude.  The  presence  of  the  nearly  free  electron  distrib¬ 
ution  in  the  solid  ensures  this  to  be  the  case  and  so  to 
reduce  the  number  of  parameters  we  have  neglected  0^  from 
the  expression  for  the  total  energy  given  in  equation  (5.3). 
Further,  since  0^  decays  with  distance,  we  only  work  to  the 
first  contributing  set  of  parameters. 

i  • ->m  equations  (5.5)  to  (5.10)  we  find  that  the  first 
contributing  interactions  to  the  phonon  frequencies  arises 
from  that  of  the  bond  electron  arising  from  next  nearest 
neighbour  bonding,  which  we  have  found  to  be  greater  than 
nearest  neighbour  bonding. 


Hence,  defining 


K„ 


d0j(r) 

dr 


=  a//2 


and 


K2  =  die2 <r) 
dr2 


=  a//2 


(5.12) 


(5.13) 


where  a  is  the  cubic  spacing,  we  find  that  the  elastic  cons 
tants  are  given  by  the  equations 
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The  Cauchy  discrepancy  is  then 

c  -  c 

12  44 

=  -  K2/(8a)  +  e” 

f 

+  2e 

V 

(5.18) 

We  have  calculated  the  phonon  frequencies  of  Fe  using 
equations  (5.10)  and  (2.29)  in  the  three  main  symmetry 
directions.  The  parameters  and  were  varied  until  the 
experimentally  determined  dispersion  curves  were  well  rep¬ 
roduced.  The  parameter  d  in  the  volume  dependent  energy  was 
chosen  so  that  the  Cauchy  discrepancy  remained  positive.  In 
figure  5.3  (a)  we  have  plotted  the  dispersion  curve:;  for  the 
best  values  of  K^,  K2  and  d,  and  in  figure  5.3  (b)  the  dis¬ 
persion  curves  calculated  with  K^,  and  K2  equal  to  zero  are 
shown.  '  comparison  of  the  two  figures  shows  the  effect  of 
introducing  the  bond  charges.  From  figure  5. 3. (a)  we  see 
the  agreement  with  experiment  is  very  good  in  the  (110) 
direction  and  for  small  k  in  the  (100)  and  (111)  directions. 

However,  near  the  zone  boundary  in  the  (100)  and  (111) 
directions,  the  agreement  is  less  satisfactory .  Also  we 
find  a  crossing  of  the  longitudinal  and  transverse  branches 
in  the  (100)  direction,  which  is  not  observed  experimentally. 
The  crossing  exists  even  with  central  forces  alone  and  hence 
it  appears  to  be  an  artifact  of  the  central  pair  potential 
that  we  have  used.  A  variation  of  the  parameters  does  rem¬ 
ove  the  crossing,  but  at  the  expense  of  agreement  in  the 
other  directions. 

We  have  also  calculated  the  second  order  elastic  cons- 
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tants  with  and  without  the  angular  contributions.  From 
equation  (5.18)  it  is  clear  that  the  angular  interaction 
gives  a  negative  contribution  to  the  Cauchy  discrepancy. 
The  values  of  the  parameters  together  with  the  elastic 
constants  are  tabulated  in  table  5.1. 


§5.4  Discussion 


As  mentioned  in  Chapter  3,  unlike  in  the  f.c.c.  metals, 
the  experimentally  determined  phonon  dispersion  curves  of 
the  b.c.c.  transition  metals  bear  little  resemblance  to  those 
of  the  b.c.c.  simple  metals.  From  our  calculation  of  the 
phonon  dispersion  curves  of  Fe,  we  find  that  positive  values 
of  and  lead  to  an  increase  in  all  the  w(k)  at  the  zone 
boundaries  and  at  the  trough  region  of  the  [ill]  longitud¬ 
inal  branch.  In  contrast,  the  effect  of  Ey(V)  is  to  increase 
the  frequencies  of  the  longitudinal  branches  whilst  decreas¬ 
ing  those  of  the  transverse  modes.  It  also  deepens  the  trough 
in  the  Qll]  longitudinal  mode.  Hence,  the  introduction  of 
the  bond  charges  provides  an  explanation  of  the  different 
phonon  dispersion  relations  observed  in  the  b.c.c.  transit¬ 
ion  metals.  Chen  and  Brockhouse  (1964),  Shaw  and  Muhlestein 
(1971),  Brockhouse  et  al  (1967),  Powell  et  al  (1968)  and 
Collela  and  Batterman  (1970). 

From  equations  (5.15)  to  (5.18),  we  see  that  the  incl¬ 
usion  of  the  angular  forces  also  provides  the  possibility 
of  a  negative  Cauchy  discrepancy.  A  positive  or  negative 
value  of  1<2 ,  from  equation  (5.18)  decreases  or  enhances 
respectively  the  Cauchy  discrepancy.  Hence,  from  the  exper¬ 
imentally  observed  e]  'tic  constants,  it  would  appear  that 
in  Fe,  and  to  a  greater  extent  in  Cr,  Kj  is  positive  whereas, 
it  is  negative  in  the  other  b.c.c.  transition  metals.  This  is 
consistent  with  our  model.  Since  Kj  is  the  second  derivat¬ 
ive  of  the  bond  electron-ion  potential,  that  is  in  all  prob¬ 
ability  oscillatory,  at  distance  a//2,  its  values  and  sign 
depends  on  the  lattice  spacing,  a.  Thus,  as  Fe  and  Cr  have 
lattice  spacings  of  less  than  3$  and  the  other  b.c.c. 
transition  metals  have  lattice  spacings  greater  than  3&,  a 
change  in  sign  in  the  second  derivative  is  clearly  possible 
if  the  wavelength  of  the  oscillatory  potential  is  of  the 
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TABLE  5.1 


The  parameters  used  in  the  calculation  of  the  phonon 
frequencies  of  Fe  and  the  second  order  elastic  constants. 

K1  =  0.1875  eV/A 
K2  =  0.0625  eV/82 

Volume  dependent  Energy  density  =  0.036  eV/X^ 


Cij(1010N/m2) 

t  n-  (a) 

Expt . 

Ev(V)  = 

k^=k2=o 

k1=k2=o 

with 

angularity 

C11 

23.3 

27.85 

31.16 

31.34 

C12 

13.6 

13.19 

14.79 

14.75 

C4  4 

11.8 

13.19 

14.15 

14.15 

(a)  Gschneider  (1964). 
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order  of  that  of  the  central  ion-ion  pair  potential.  Hence, 
this  model  provides  the  facility  of  calculating  a  negative 
Cauchy  discrepancy  in  Cr. 

From  our  calculation  of  Fe,  it  is  clear  that  a  good 
pair  potential  is  required  before  the  effects  of  the  angul¬ 
arity  in  the  charge  distribution  can  be  studied.  Figure 
5.1  shows  that  the  central  ion-ion  pair  potential  that  we 
have  extracted  from  the  liquid  state ,  although  fairly  good 
in  reproducing  the  experimentally  determined  phonon  disper¬ 
sion  curves  of  Fe,  introduces  a  crossing  of  the  longitudinal 
and  transverse  modes  in  the  [jLOCQ  direction.  Although  this 
crossing  is  shifted  closer  to  the  zone  boundary  on  the  intro¬ 
duction  of  the  bond  charges,  its  persistence  leads  us  to 
believe  that  the  central  ion-ion  pair  potential  that  we  have 
constructed  may  have  to  be  further  refined.  However,  it  is 
interesting  to  note  that  even  in  the  simple  metals  Na  (Woods 
et  al.  1962)  and  K  (Cowley  and  Woods  1966)  in  which  the 
force  fields  are  central,  this  crossing  is  observed  and  thus 
we  can  safely  assume  that  it  is  not  a  feature  of  the  angul¬ 
arity  in  the  force  fields. 

W"  ^ive  thus  shown  that  the  bond  charge  model  described 
in  §5.1  is  very  useful  in  understanding  the  phonon  disper¬ 
sion  relations  of  the  transition  metals.  Further,  .the  gen¬ 
eral  agreement  with  experiment  leads  us  to  conclude  that 
next  nearest  neighbour  bonding  is  more  important  than  that 
of  nearest-  neighbour  bonding  in  the  b.c.c.  transition  metals. 
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CHAPTER  6 


Conclusion 


The  problem  of  determining  the  force  field  in  a  tran¬ 
sition  metal  is  clearly  bound  up  with  that  of  calculating 
the  electron  distribution.  For  example,  in  lattice  dynamics 
within  the  adiabatic  and  harmonic  approximations,  the  entire 
problem  of  the  dynamical  matrix  and  hence,  the  phonon  disper¬ 
sion  relations  hinges  on  calculating  the  change  in  the  elec¬ 
tron  distribution  as  the  ions  are  moved  small  distances  from 
equilibrium. 

Pair  potential  theory  assumes  that  the  localised  charge 
distributions  move  rigidly  with  the  ions  as  they  vibrate. 
Though  pair  potentials  afford  the  starting  point  of  all  the 
work  prec:ented  here,  the  directional  character  of 
the  d  electrons  in  the  transition  metals  demands  corrections 
to  this  description  of  the  force  fields. 

Of  course,  it  is  important  in  the  d-ba^d  metals  to 

recogr that  the  electrons  in  these  materi  als  have  widely 

varying  properties  as  one  goes  from  itinerant  s  electrons 

to  spat nally  localised  d  electrons,  with,  naturally  enough, 

hybrid s-d  electrons  also  playing  an  imp-rtant  role  in 

the  charge  density  and  hence  the  force  field'’. 

•  • 

Thus,  incorporating  the  unhybridized  d  electrons  into 
the  ionic  charge  distribution,  the  picture  of  a  transition 
metal  is  that  of  a  collection  of  ions  situated  at  lattice 
sites  and  localised  electron  distributions  arising  from 
s-d  hybrid  overlaps  immersed  in  a  nearly  free  electron  gas. 

In  Chapter  2  we  saw  that  the  total  energy  of  a  crystal  could 
be  written  as  a  sum  of  ion-ion  pair  interactions.  In  the 
above  picture,  the  ions  thus  interact  with  each  other  direc¬ 
tly,  via  the  electron  gas  and  via  the  localised  electron 
distributions.  We  also  saw  that  it  is  the  last  that  gives 
rise  to  angularity  in  the  ion-ion  interactions. 


To  answer  the  first  of  the  two  questions  posed  in  the 
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Introduction,  we  see  that  a  central  pair  potential  descrip¬ 
tion  would  be  valid  only  in  those  metals  where  the  local¬ 
ised  electron  distributions  are  very  small.  From  the  argu¬ 
ments  presented  in  Chapter  2  we  saw  that  this  would  be  best 
for  those  transition  metals  possessing  the  f.c.c.  structure 
The  usefulness  of  the  central  pair  potential  for  Ni, 
extracted  from  the  liquid  structure  data,  in  reproducing 
the  experimentally  determined  phonon  dispersion  curves, 
suggests  a  starting  point  in  the  calculation  of  a  central 
pair  potential  in  solids.  This  view  is  confirmed  by  our 
calculations  on  Pt.  The  liquid  structure  calculated  from 
the  central  pair  potential,  obtained  by  fitting  the  phonon 
dispersion  curves,  shows  very  good  agreement  with  exper¬ 
iment.  Thus,  in  Chapter  2  we  showed  that  the  concept  of 
central  forces  is  a  useful  one  in  the  f.c.c.  transition 
metals  and  the  central  pair  potential  in  solids  is  directly 
related  to  that  in  the  liquid  state. 

In  an  effort  to  obtain  information  about  the  angular¬ 
ity  in  the  charge  distribution  and  hence  the  force  fields, 
we  investigated  the  available  X-ray  scattering  factor  data 
on  the  transition  metals  in  Chapter  3.  By  formulating  the 
theory  in  terms  of  WF's.  we  showed  that  the  charge  distrib¬ 
ution  in  the  transition  metals  could  be  represented  by  bond 
charges  situated  mid-way  between  nearest  and  next  nearest 
neighbours.  The  validity  and  usefulness  of  this  formalism 
was  demonstrated  with  a  calculation  on  Be  in  Chapter  4. 

Thus,  with  a  useful  representation  of  the  charge  den¬ 
sity  in  the  d-band  metals  we  proceeded  to  construct  a  model 
by  which  the  phonon  frequencies  of  Fe  were  calculated.  In 
this  model,  presented  in  Chapter  5,  the  force  fields  consis¬ 
ted  of  bond  electron-bond  electron  and  bond  electron-ion 
interactions  in  addition  to  the  effective  ion-ion  interaction 
present  in  the  simple  metals.  Confidence  in  this  approach 
has  resulted  from  the  application  to  the  calculation  of 
the  phonon  dispersion  relations  in  Fe,  where  fairly  satis¬ 
factory  results  were  obtained.  The  ion-ion  interaction  in 
the  transition  metals  has  the  same  form  as  in  the  simple 
metals,  i.e.  given  by  equation  (2.8).  From  our  work  we 
see  that  a  semi-empirical  way  in  which  we  can  determine  the 


parameters  is  by  appealing  to  the  liquid  structure  factor 
data.  The  repulsivity,  n,  can  be  obtained  from  an  analysis 
of  melting  and  radiation  damage  data. 

The  volume  dependent  energy,  Ev(V)  and  the  angular 
contribution  can  be  represented  in  some  form  and  the  para¬ 
meters  determined  by  relating  it  closely  to  experimental 
data. 

By  perturbation  theory,  the  bond  electron  interactions 
will  also  be  oscillatory  in  nature.  However,  it  is  not 
clear  whether  it  is  correct  to  consider  all  the  two  body 
interactions  separately.  A  more  correct  approach  would  be 
to  treat  all  the  electrons  together  in  a  mixed  WF  and  plane 
wave  formulation.  However,  from  this  first  principles 
approach  as  formulated  by  Pick  et  al.  (1970)  and  Hanke 
(1973),  it  is  not  easy  to  understand  the  nature  of  the 
force  fields  in  real  space. 

Finally,  we  wish  to  make  some  comments  about  the  acc¬ 
uracy  of  the  potentials  that  have  been  constructed  in  this 
thesis.  Though,  as  mentioned  in  §2.6  the  parameters  of  the 
central  pair  interaction  are  accurate  to  better  that  95%, 
it  must  be  noted  that  we  have  assumed  a  particular  form  for 
the  interaction.  This  form  was  motivated  by  that  arising 
fr'"-  the  asymptotic  form  of  the  displaced  charge  around  an 
ion  Ln  any  electron  gas,  and  is  rigorously  valid  only  at  very 
large  distances.  We  have  in  this  work  also  not  taken  into 
account  the  changes  in  the  potential  due  to  density  and 
pressure  variations,  although  we  recognised  that  the  cen¬ 
tral  pair  potential  in  the  liquid  differs  from  that  in  the 
solid.  From  figure  2.2.,  we  see  that  again  the  differences 
are  only  of  the  order  of  5%.  The  aim  of  this  work  has, 
however,  been  to  obtain  practically  useful  representations 
of  the  force  fields  in  the  transition  metals.  The  results 
which  have  been  obtained  by  appealing  to  certain  experimental 
data  must  be  able  to  reproduce  other  quite  different,  exper¬ 
imental  results  before  any  claims  of  a  true  representation 
of  the  force  fields  in  the  transition  metals  can  be  made. 

For  example,  an  important  question  for  future  work  would 
be  whether  or  not  such  force  fields  were  useful  in  calcul- 


ating  vacancy  and/or  self  interstitial  energies  in  the 
transition  metals.  This  problem,  and  more  generally  that 
of  defect  interactions,  is  of  sufficient  technological 
importance  to  warrant  a  fully  study  with  the  currently 
available  force  fields. 
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